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ABSTRACT 



The equations for estimating the performance of an axial flow 
turbine are examined to determine the efficiency and other parameters 
when very high power outputs per stage are required. It Is Indicated 
that turbine stages should be designed for peak efficiency or for 
high power output, but that compromise designs are not advisable. 

A survey of the Immediately available literature on losses In 
turbine stages at the U. S. Naval Postgraduate School Is outlined. 

An attempt to find and correlate specific loss data with theoretical 
relations was not successful. Data Is fragmentary, theoretical re- 
lations are few, and the correlation of the available data on profile 
losses, tip clearance losses, effect of trailing edge thickness and 
secondary losses Is poor. 

An axial flow single stage turbine was designed for transonic 
flows and a digital computer program was written to evaluate the 
off-deslgn performance of this turbine. The particular relations for 
supersonic expansion after a blade row are Included In the program 
as well as the effects of entry shocks Into the rotor. Initial 
Indications are that the program should be useful for evaluating 
turbines designed by the procedure shown. 

The computer programs Included are written in FORTRAN language 
for the CDC 1604 Digital Computer. 
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THE PERFORMANCE OF HIGH POWER OUTPUT AXIAL 



FLOW TURBINES UTILIZING FLOWS IN THE TRANSONIC REGIME 

I. Introduction 

The increasing need in space applications for turbines of high 
power output from a minimum physical volume and weight requires more 
knowledge of efficiency, losses and performance estimation of these 
turbines • 

If the required power output and space limitations are predominant 
then efficiency must suffer. As a means of showing the degree of 
efficiency sacrificed, the equations which determine the power output 
and efficiency of an arbitrary turbine were programmed for digital 
computer solution in terms of the head coefficient. The head coef- 
ficient for peak efficiency was determined and then the head coef- 
ficient was increased in increments and various parameters affecting 
the turbine performance were computed at each increment until the 
efficiency dropped to an arbitrary cutoff point. The data are plotted 
in non-dimensional form and simple equations are found to represent 
the data. 

Of paramount importance in the design of a turbine is the esti- 
mation of the losses due to boundary layer formation, mixing effects, 
clearances and secondary flows. This area is explored in the hope 
of finding a more accurate and precise method of presenting these 
losses. A survey of the immediately available literature is included. 



A particular single-stage turbine is designed for supersonic 
flow leaving the stator blade row and a relative subsonic flow as 
seen by the rotor. The equations for determining the effects of 
after-expansion behind a blade row and entry shocks entering a 
blade row are programmed for the computer in addition to the 
regular performance parameters. The rpm and pressure ratio across 
the turbine are then varied to show the effects of demanding extreme 
power outputs from a turbine so designed, A performance map of the 
turbine is developed. 
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II. Development of High Head Coefficient Data 



A. Basic Equations 

In Ref. 1 the equations are developed to establish the design 
parameters for turbine stages. These equations are obtained by 
assuming that the flow through the turbine is represented by the 
conditions on a mean stream surface and the primary parameter is 
taken to be the "head coefficient". The head coefficient relates 
the total isentropic drop of the stage from the inlet stagnation 
pressure to the discharge static pressure to the kinetic "energy" 
of the peripheral speed at the rotor entrance and is defined as 




( 1 ) 



In order to relate the theoretical expressions to values found 
in engineering use, some expression for the losses involved must be 
used. A velocity coefficient ( ^ ) is used as a measure of rotor 
efficiency and is defined as 




( 2 ) 



This velocity coefficient has been found through experience to depend 
primarily on the deflection of the flow in the rotor ( ) . Vavra 
has established a mean curve through data from a variety of sources 
to be 






4S7 



( 3 ) 
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and this equation is used to represent the losses for purposes of 
evaluating the efficiency at each step of the calculation. 



Following the procedure outlined in Ref. 1 , arbitrary values 
were chosen for the angle of the flow leaving the stator y 

reheat factor (f), and the carryover coefficient of kinetic energy 
into the next stage the radius ratio (R2/R^) , the meridional 

velocity ratio ^ velocity coefficient in the stator (fy. 



and the carryover coefficient of kinetic energy from the stator into 

* 

, . TVipn Hpat-pp ^“ppr^^nn 

the head coefficient is increased in increments. The degree of 



the rotor Then the degree of reaction (r ) is selected and 



reaction is the fraction of the isentropic enthalpy change through 
the turbine stage which occurs in the rotor. At each increment of 
head coefficient the particular values of velocity coefficient y 

internal efficiency ( > leaving loss coefficient (kg) > ratio of 
meridional absolute velocity to peripheral speed y relative 

velocity ratio (W2/Wj^), and the relative and absolute angles (0(^2> 
/^l’ 1 ^ 2 ^ computed. 



A similar computer program using the same equations was written 
to find the peak efficiency for a given entering angle and degree of 
reaction, and then all the other parameters are computed. The values 
are also found for higher head coefficients such that the efficiency 
is 0.57o less than optimum and then 1.07o less than optimum. 
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B. 



Computer Program 

The equations to be solved by iteration until are; 



1 . 



• 



/ 









4' 



X 



where 



(^z/^/) /kf(LO(^ ''J 
/^r= (/+f)yX^+ 

-2(4 yWflT, 



( 4 ) 



(5) 



(5a) 



3- = 



4.97 



/o^ r /8o 



(3) 



Once this condition ( ^ * ^ ) is satisfied the values of the 

A 

other parameters are found from 



1. 



•J 

4 = -4 



( 6 ) 



(7) 




'^\ , / 






c 

U/ 



= f -jfh^ 'ffX ^o^o(i 



( 8 ) 
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■A / 'f^Ai 



(9) 




The equations were coded in Fortran language, Ref. 2, and solved 
on the Control Data Corporation 1604 computer at the U. S. Naval 
Postgraduate School. A copy of the program is presented in Appendix III. 
The resulting data printouts are included in Appendix V as Table I. 



The program finding the peak efficiencies, efficiencies 0.5% less 
than optimum and efficiencies 1.07« less than optimum is also in 
Appendix III, with data in Appendix V as Table II. It should be 
noted here that all equations are completely coded in Fortran language 
so that any selected numerical values may be introduced in the be- 
ginning of the program. Any of these arbitrarily chosen values may 
be changed and the equations of the main program remain unchanged. 



C. Results 

Representative data from Tables I and II are presented in Fig. 1. 
The regularity and progression of the data curves indicate that non- 
diraensionalizing the efficiency and head coefficient parameters might 
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I 



t'S profitable. When all the efficiencies for one entering angle and 

one degree of reaction are divided by the optimum efficiency for that 

set of conditions (^/ the head coefficients are divided 

by the optimum head coefficient (k. /k. \ the curves of Figs. 2. 

IS IS opt ° ’ 

3, 4 and 5 result. 



A closer inspection of these plots reveals that for each degree 
of reaction, one curve suffices for all four entering angles. 

The cubic equation 




-.000223 



/‘"is ^ 


3 

+ .008345 1 


I k ^ 

^ is 


2 

T O ^ O i 


Ihs \ 


^ is opt 


k 

i is opt 


- . 12068 


1 k. 

1 IS opt j 



+1.125 



represents all the data for r “0.0 for all entering angles from 60° 
to 75°. 



Similarly as an approximation to all the data at r “ 0.25, the 
equation is 






k. 



is 



is opt 



+.00770 



k. 

IS 



k. 



1098 



k. 

IS 



IS opt 



k. 

IS opt 



+ 1.1 



and for r * 0.50 




-.000125 



/k. 

IS 



k. 

IS opt 



+.00554 



is 



k. 

IS opt 



-.09414 



k, 

IS 



k. 

IS opt 



+1.097 



The maximum error in the range of entering angles presented here is 
approximately 8 parts in 130. This is probably accurate enough for 
preliminary engineering estimates based on the equations in Ref. 1 
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used in the computer programs. As will be seen later, the loss data 
is not known to this accuracy. 

The curves show that, roughly speaking, increasing the head 
coefficient to five times the optimum causes a reduction of 307o of 
the obtainable peak efficiency. Increasing the head coefficient to 
ten times the optimum causes a reduction of 457o of the obtainable peak 
efficiency. This indicates that if it is necessary to go to very high 
head coefficients to obtain a certain required work output, the rate 
of change of efficiency with increasing head coefficient is greatest 
from optimum to five times optimum head coefficient. The rate of 
change of efficiency decreases as the head coefficient is further 
increased. 

Table I of Appendix V can be used to establish the other stage 
parameters after the head coefficient range of interest is chosen. 

Figs. 1 through 5 are not the off-design performance of a 
particular turbine, but they are the locus of points of an infinite 
number of different turbine designs. Each point on each curve 
represents a turbine designed for that condition. Therefore, if the 
requirement is that the turbine be designed for a high head coef- 
ficient, the loss in efficiency need not be the absolute controlling 
factor, and the design may as well be based on extreme high values 
of head coefficient. These investigations show that it is advisable 
to either design the turbine for peak efficiency of maximum output 
but do not compromise both by operating at only two or three times 
the optimum head coefficient. 



- 8 - 









J J 



fm 



III. Loss Investigation 



Since the design of a turbine depends on the loss assumptions 
made, the use of an average velocity coefficient based on general 
considerations of flow deflection in the rotor is not a very precise 
method of estimating the losses. This was the method used in the 
preceding section of this thesis. A better way should be found and 
the following considerations are basic to the problem. 

A. Loss Components 

It seems logical to class the losses encountered in turbines 
into four particular categories; 

1) Profile loss -- the friction loss associated with the 
formation of the boundary layer on the blades of the cascade. 

2) Tip clearance loss -- the pressure and energy loss associated 
vath the leakage flow between the tips of the blades and the casing 

of the rotor, including the effects of trailing vortices at the hub 
where there is zero clearance. 

3) Trailing edge loss the loss due to mixing effects behind 
the blades of the row as caused by the finite thickness of the blade 
trailing edge. 

4) Secondary loss -- the loss due to the vortices which are set 
up in the flow in the curved channel between blades giving a velocity 
component normal to the desired direction. 
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In the axial flow turbine, the closely packed blades and large 
flow deflections usually require the flow through the turbine to be 
considered from the point of view of ''channel flow" between blades 
instead of using airfoil theories. As will be seen in the subsequent 
discussion, this secondary flow causes the principal problem pre- 
venting the accurate assessment of the losses. 

Looking at the flow between two blades of an arbitrary cascade. 
Fig. 6, the velocity distribution is as shown for the mid-span of 
the blade channel and it can be seen that a static pressure gradient 
is developed between surfaces 0-1-2 and 0*-l*-2*. The centrifugal 
forces on the fluid being turned cause an increase in the static 
pressure near the 0*-l*-2* wall and decrease the static pressure 
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Assuming non-viscous flow outside the boundary layer requires the 
flow to increase velocity near the 0-1-2 convex wall to maintain the 



velocity must decrease near the 0*-l*-2* concave wall outside the 
boundary layer. In the boundary layer at the top and bottom of the 
channel, the loss of velocity (and total pressure) due to friction 
causes the flow direction to be dominated by the pressure gradient 
imposed (assuming the static pressure is constant through the boundary 
layer). Therefore, as the flow progresses from station 0-0’ to 
station 1-1* to station 2-2* , the fluid in the upper and lower 
boundary layers tends to flow from 0* toward 1 and from 1* toward 2. 

This change of direction of the flow in the boundary layers 
results in the formation of two vortices superimposed on the main 
flow leaving the channel at station 2-2*. In Fig. 7 the flow patterns 
in the upper and lower boundary layers are shown as well as the 
induced vortices rotation direction at the exit plane. 

This vorticity results in velocity components developed perpen- 
dicular to the desired exit velocity. Since a velocity component 
normal to the desired direction cannot be recovered usually, most of 
the useful energy in these components is lost in friction. This is 
the so-called **secondary loss**. For a rotating row of blades, the 
blade tip clearance also has a marked effect on the secondary flow. 

The leakage from the high pressure to the low pressure side of a blade 
encourages an additional trailing vortex to form a reaction blading. 



total 




Conversely, the flow 
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For impulse blading the flow over the tips which is not deflected also 
enhances vortex fonsation. For cascades of infinite blade height 
(two-dimensional only) there is no secondary flow since the mechanism 
depends principally on the annulus boundary layers. 




With the physical parameters of blade height, trailing edge 
thickness, tip clearance, profile shape, and spacing given, it should 
be possible to make a reasonably exact estimate of losses. 

B. Survey of the Literature 

ft 

The iaDediately available literature at the U. S. Naval Post- 
graduate Sphool was surveyed to find the particular loss breakdown 



used by the various sources, the theory and representations offered, 
and any precise loss data published. The survey is limited in scope 
and critical since the particular desire is to locate specific loss 
theories and data. The results of this survey are presented here 
in the order of reference books, general texts, and then periodicals 
and reports. The NACA/NASA reports provide a wealth of data on the 
efficiency, or total loss coefficient, of specific turbine designs 
from tests. Very little information was found for loss components, 
however, in the frame of reference considered here. 



1. In Ref. 1 an equation is derived for the determination of 
a flow coefficient in a rotor as a function of the change in 

flow angle of the relative flow. This equation is a fitted curve 
through data from several sources and a variety of blade shapes and 
operating conditions and is the same equation used in the preceding 
section of this thesis. In Fig. 8 is shown the general shape of the 
curve and the derived equation (3). 



A loss coefficient can be defined as 



Considering the case of 4^ =0, we see the loss coefficient 
’C ■ ‘ - /'o represents those losses in the cascade which are 
present when the flow is not deflected. As the flow is deflected 
an amount y the reduction in ^ indicates an additional loss. 

The additional loss will be taken as the secondary loss. The 
secondary loss coefficient is then 






( 11 ) 



( 12 ) 
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Curves shown in subsequent discussions of secondary losses will now 
show a loss coefficient curve based on this equation plotted as a 
dashed line and labeled ”VAVRA**. This curve will act as a general 
reference and will show **average” conditions so a comparison of each 
method can be made. 




Fig. 8 Flow Coefficient as a Function 
of Turning Angle in a Cascade 

2. A veritable fountainhead of theory, design and ideas are 
presented in Ref. 3. This basic work by Stodola touches on almost 
all the subjects considered important now. The important works of 
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others up until 1927 are included to show their measurements and de- 
velopments. The losses through fixed rows of blades are given in 
the form of the so-called velocity coefficient ^ ; the ratio of actual 
velocity to theoretical velocity. The measurements of the Institution 
of Mechanical Engineers of England and others all show good correlation 
for flows up to sonic speed for stators with an axial entering angle 
and an exiting angle of 70°-78° from the axial. For supersonic 
speeds the data show some scatter due to difference in measuring 
apparatus and some basic differences in design of the convergent- 
divergent passage. 

Secondary flows in bends are explained but there is wide vari- 
ation in the available data for friction losses in bends. No pub- 
lished data were available for the losses in rotating blade rows. 

The data and experiments on fixed rows of impulse rotor blades were 
of limited application and only show general trends. 

For design purposes, velocity coefficients and loss coefficients 

y 

are given in general terms only, i.e., ‘‘velocity coefficient for 

I 

stators * 0.975-0.92 for long and short nozzles respectively“ , The 
curve presented by Stodola for the velocity coefficient for rotating , 
blade rows as a function of turning angle is included in the data 
used for the determination of the equation (3) by Vavra . 

The tip clearance losses had been explored in a limited fashion. 
These losses were accounted for by an empirical formula for the de- 
crease in internal work as a function of clearance and blade height. 
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( 13 ) 



- 6.27 k^'^/h 

where k * radial clearance and h is the blade height corresponding 
to zero clearance. This relation, credited to Anderhub, is 
plotted on Fig. 9 and noted as . 

The steam turbine design procedures presented are strongly 
dependent on making the new design a small modification of an already 
designed unit which is in service and for which the efficiency 
is known. The peak efficiencies were recognized to be in the regions 
of k^^ from 1.0 to 5.0, as shown in the previous section of this 
paper. 

3. The axial turbine section of Ref. 4 gives a lucid and brief 
explanation of the flows in such machines but the loss data and 
secondary flow considerations are based completely on Ref. 22, 23 
and 24 which are considered later in this paper. 

4. Since books written to be used as school texts of necessity 
try to cover the complete spectrum of fluid mechanics, theory, and 
design, only the most cursory reference is made to specific losses 
in turbine stages. Refs. 5, 6 and 7 are representative of books of 
this type. Ref. 5 defines a work recovery coefficient as the ratio 
of useful work extracted by the rotor to the total kinetic energy 

of the flow leaving the stator, times the stator flow efficiency. 

This work-recovery coefficient variation with the velocity ratio 
(U/V^) is shown for impulse and reaction conditions and no other 
explanation of the nature of the losses is offered. The charts 
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presented are based on data from Ref. 16, in which report the losses 
in a blade row are presented only as a function of stagger angle and 
incidence. 

In Ref. 6 the author makes the statement that no adequate 
criterion is available for the difference in the flow conditions 
between an efficient and an inefficient hydrodynamic machine of a 
given type. Any losses and efficiencies of turbines offered are 
tied to airfoil theory or definitions of losses by Stodola and 
others. Secondary flows and radial flow of the boundary layer are 
mentioned only from Ref. 17 where the measurements were made on 
widely spaced blades on a rotating blade row at low speed. 

A large list of references on the work in gas turbines is 
given in Ref. 7 for the data available at that time. Principally 
the work in Ref. 16 below is mentioned in connection with turbines 
but only general terms are offered and no specific loss derivations 
are offered. 

5. The most comprehensive effort to make a breakdown of losses 
and show specific data found in a textbook is that of Ref. 8. This 
book takes selected data from Refs. 15, 16 and 23, discussed later 
in this paper, and presents plots of the losses due to leakage, 
profile losses, losses from turbulence and wall friction (including 
secondary loss effects) and incidence losses. The discussion is 
brief, however, and no new data or theory are presented. 
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6. An attempt to mathematically tackle the problem of the 
secondary flows in bends is presented in Ref. 9. The equations of 
motion of an ideal fluid are used and frictional effects are ignored. 

An induced drag pressure drop coefficient is defined to account for 
the power loss due to secondary flow. The method is applied to the 
flow in the corner of a 3* x 4* wind tunnel which has two turning 
vanes to simulate a blade row and the method shows the proper trends. 

A quantitative value derived for the loss is computed and found to 

be about twice the measured loss. No other example of the theory 
being applied is available. 

7. A qualitative investigation of the flow in boundary layers 
and wakes of blade rows is offered in Ref. 10. An excellent discussion 
of the general nature of these flows is accompanied by the statement 
that a quantitative theory for these effects is not available. 

8. An analytical method of estimating turbine performance 
described in Ref. 11 is based on losses due to incidence angle and 
losses due to flow of the fluid through the blade row. A blading 
loss parameter is used to represent the losses but only a general 
range of values of the parameter are given. No method is presented 
for determining the magnitude of the blading loss parameter before 
the turbine is built. 

9. The general state of the art to 1948 is described in Ref. 12. 
The paper is not intended to be used as a basis for other than general 
design considerations, however, and some empirical relations presented 
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were attempts to cover both axial turbines and compressors. This 
necessarily limited the worth of such relations for large flow 
deflection turbine loss estimations. Those relations presented for 
strictly turbine use are repeated in Refs. 24 and 25 which will be 
covered in more detail later in this thesis. 

10. The losses in stators are, of course, easier to measure 
than the losses in rotors. Some excellent tests and measurements 
are presented in Refs. 13 and 14. These show that, in general, a 
properly designed nozzle for a turbine will have a velocity coef- 
ficient (ratio of actual velocity to isentropic velocity) value of .96. 
This single value is found to be uniform over a considerable range 

of designs and conditions. Off-design values do fluctuate, partic- 
ularly in the case of supersonic nozzles but the reasons do not seem 
obscure and loss estimations for stator rows should be on relatively 
firm ground. 

11. Ref. 15 is mentioned as a reference work in many textbooks 
on turbomachinery. This compilation of the then current design 
practices only shows the losses as a function of axial length of the 
blades of the rotor, turning angle in the rotor and incidence angle. 

The curves for loss as a function of turning angle are included in 
the data used for the determination of the equation (3) by Vavra. 
Specific detailed loss breakdowns or data are not available. 

12. The systematic investigations into secondary loss presented 
in Refs. 18 and 19 are of real interest although the test data is too 
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limited to permit any universal application. An attempt is made to 
define the two-dimensional loss by theory, measure the end loss at 
the blade root and the clearance loss at the tip, and subtract these 
from the total pressure loss measured to arrive at the secondary 
loss. The tests are performed on a particular airfoil shape in 
three solidity arrangements and three flow deflections. The inac- 
curacies and errors discussed should help any further investigator 
in this subject. 

13. A comparable breakdown of the losses into components as 
assumed here is offered in Ref. 20. This work by Markov also includes 
the experimental results of some other Russians in the field. With- 
out attempting to completely describe the development of equations 
and data presented, the pertinent data for the present consideration 
is extracted. 

The profile loss, or two-dimensional loss which would be caused 
by blades of infinite length is shown on Fig. 10 along with other 
data to be described later. The basis of the Ref. 20 data is a whole 
series of experimental measurements of widely used cascades for 
turbines in Russia. For stators, the profile loss coefficient is 
almost constant at .03 over a wide range of conditions. 

The tip clearance loss is shown on Fig. 9 as represented by two 
plots in Ref. 20. One plot is shown as the loss for a reaction stage 
and one plot is based on experimental data for the efficiency of stages 
with a small degree of reaction. The magnitudes of the degree of 
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reaction are not specified but the trend is correct. Stages with a 
high degree of reaction have a higher loss across the tip, as would 
be expected due to the pressure difference across the rotor from 



entrance to exit. 

The trailing edge loss for cases where the pressure gradient 
is small along the surface of the trailing edges of the blades is 
computed by 



Uhere such a pressure gradient is present the use of boundary layer 
thickness parameters is required. The trailing edge loss coefficient 
according to eq. (14) is shown in Fig. 20. 

The energy lost due to the secondary flow is treated as a 
function of the mass flow rate in the cascade. The efficiency of a 
cascade of blades with finite length is defined as 




(14) 




(15) 



The loss coefficient is defined as 




(16) 



and so 
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The secondary flow loss coefficient 
used for comparison in this thesis by 



is converted to the 




^ 5^ ^(2-dim) 

An expression for ^ is developed that is a function of the circu- 
lation around the blades, the blade width, and incompressible flow, 
namely 







where the coefficient K depends on the expansion of the flow in the 
cascade. The value of K is determined from experimental data to be 



K ^0.0/ 






Thus 



f. = m 



L J 



(17) 



In Fig. 11 the function ^ is shown for the two-dimensional efficiency 
of Ref. 1, 

^ =■ ^^ = 0.9Z6Z 

/(2-dim) fo 



Values of h/c plotted are chosen arbitrarily but Ref. 20 specifies 
that the derivation holds for blades that are not excessively short. 
There is assumed to be some two-dimensional flow in the channel be- 
tween the blades at the center section. Experimental data are shown 
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for short blades, and these data commence at 




■ 2.0 and go to 



smaller values. So the plot in Fig. 11 is based on eq. (17) and 



resentative values based on the Vavra equation (3), but these are 
specific data for short blades not arrived at using the equation 
above. The secondary loss as computed by Markov is low compared 
to any other available data. When complete loss coefficients are 
computed in Ref. 20 the values are comparable to any other method 
however, since the tip clearance losses, trailing edge losses, 
blade height corrections, etc., used bring the total loss to the 
proper level. Any use of the Ref. 20 data must take account of 
this fact. 

Of interest also in Ref. 20 is a relation for the discharge 
flow angle of a turbine blade row. The relation is 



and it gives more accurate results than the other widely-used 
equation 



is shown for a value of |*“| ■ 2.0, the lower limit permitted. For 
[ ■ 1.33 and 1.25 the mean coefficients are about 107o of the rep- 



c I 




(18) 




(19) 



which does not take account of trailing edge thickness. 
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14. The tip clearance loss as represented by Ref. 21 is also 



shown in Fig. 9. The equation derived is 







( 20 ) 



where 



- .1011 + 4.667 (k/b) 

- -.0422 + 2.790 (k/b) 

£ 



h • blade height 



b ■ blade width 



k ■ clearance 



and the results are shown for two values of (h/b). 

Also given in Ref. 21 is the equation for the minimum induced 
drag due to flow in a clearance gap at the tip of an airfoil (as 
derived by Betz) to be 



where K 2 is a form factor depending on the clearance, blade spacing 
and relative exit angle of the flow. 

Assuming the axial component of velocity to remain constant 
and drawing an arbitrary turbine stage velocity triangle (such that 
all angles are positive), Fig. 12a, the circulation ( /^ ) can be 
computed by integrating around the closed curve in Fig. 12b to give 



Z 





( 21 ) 




( 22 ) 
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Then the lift per blade is 



^ ( 23 ) 



The lift and drag forces exerted on a blade can be seen in Fig. 13 
as well as the axial and peripheral components of the resultant 
force F. It can be seen that for each blade channel 



The secondary flow loss will be associated with the loss from the 
pressure drop in the axial direction. 

To arrive at a loss coefficient expression for the axial force 
component, it is necessary to define the loss coefficient in terms 
of the pressure loss, as shown in Fig. 14. The loss coefficient is 
defined as 




(24) 




(25) 




(26) 



Substituting 
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and for incompressible flow. 






and collecting terms, 




Then from Fig. 13, 





H- 



Using the above derived expressions for L and 0 and equating 



(28) 



(29) 






(30) 



gives 

For an arbitrary value of ^2(^) * secondary loss is plotted in 

Fig. 15. The value of ^2(h} adjusted to change the general 

inclination of the family of curves shown to any desired value. That 
is, a lower value of K2(h) curves and vice versa. 
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15. In Ref. 22 an attempt la made to find a mathematical 
solution to the problem of secondary flows in particular. The 
approximate relations derived are not able to satisfy the boundary 
conditions. One conclusion drawn from the derivations is that the 
secondary flows are a direct consequence of a non-uniform approach 
velocity, and, if the lift is assumed proportional to the square 
of the local velocity, then the secondary flow is independent of 
aspect ratio. These conclusions are counter to the ideas of most 
investigators in the field. 

The tests and measurements made sho%fed odd discrepancies and 
it was felt that rotor-stator induced effects may have been the 
cause. No relation was deduced which is directly useful here. 

16. The ideas and relations presented in Ref. 23 and 24 were 
carried through and refined and elaborated in Ref. 25. All three 
references have the same authors. Therefore, Ref. 25 will be dis- 
cussed here as being most useful, and parts of Ref. 23 and 24 will 
will be extracted as necessary. 

By far the most detailed method presently available for deter- 
mining loss coefficients is presented in Ref. 25. The loss is 
broken into component parts similar to those assumed at the begining 
of this discussion and the quantitative data offered are based on 
tests of turbines. 
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The loss coefficients of Ref. 25 are 



Y 




Loss of total head pressure 

Total pressure at blade outlet-static pressure at blade outlet 

(32) 



The loss coefficients used In this thesis are 



f =/ -7 




(33) 



and the necessary conversion is 




For profile losses, two figures are presented, one for nozzle 
blades and one for impulse blades. These are reproduced here as 
Fig. 16a and Fig. 16b. An interpolation method is proposed for rotors 
with reaction and a method of extending the profile loss to conditions 
of other than zero incidence is also shown. 



A number of relations for the loss due to tip clearance are 
presented in Ref. 25 as 




^ ■ Hi) 



These relations are shown in Fig. 9. 



(Stodola, Reaction turbine) (35) 
(Meldahl) (36) 
(50% Reaction turbine) (37) 
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Expressions are also derived for tip clearance loss based on 
and C If Che ratio spacing/blade height Is small and the gas turning 
angle Is small. For our uses here, the gas turning angles are large. 
The loss relation arrived at Is , 

where B ■ 0.5 for blades with clearance and B • 0.25 for blades 
with shrouds. 



Equation (38) Is used to find the range of loss coefficients 
for the Impulse conditions of cos • 1.00 and • C(^. Shown 

in Fig. 9 Is the band of values corresponding to o ( ■ 60° to 75°. 



The data available to Alnley and MaChleson led them Co Che con- 
clusion that the effect of area ratio is more Important to secondary 
flows than Is the magnitude of the gas deflection. The area ratio 

used is 



A 2 (annulus area at outlet from blade row) (cos o( ^ 

19 »ll» ■■■I I I II I ■ . ■ ■ ■■■■II ■■■!■■ ■■■■■ 

A^ (annulus area at Inlet Co blade row) (cos ^ ^) 

where o( ^ « gas outlet angle measured from the axial and • blade 
inlet angle measured from Che axial. The hub ratio (Inner diameter/ 
outer diameter) vas also felt to exert a strong influence. This led 
to the plot of (A2/Aj) /Cl + (ID.OD)] versus a parameter ^ used in 
defining the drag coefficient due to the secondary flow as 




( 40 ) 
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This plot is reproduced here as Fig. 17. Ultimately then, a loss 
coefficient for secondary flows Is presented as 



for a wide span of Inlet and outlet angles. 

For trailing edge thickness losses, the assumption Is made that 
the previously derived losses are for a blade of trailing edge thick- 
ness which Is 2% of the blade pitch. For other trailing edge thick- 
nesses a correction is applied according to a plot which Is reproduced 
here as Fig. 19. 

C. Results of Loss Investigation 

The results of the loss Investigation will be discussed according 
to the original breakdown of loss components mentioned at the beginning 
of this section of this thesis. The specific turbine design used In 
the last section of this report will be used as the object of a 
quantitative comparison of the loss relations shown here. 

1. The profile loss estimation method presented by Ainley in 
Ref. 25 is straightforward and complete and should give good data 
for the "conventional" blade shapes from which the data is derived. 

An anomaly presents Itself In Fig. 10 however, when a comparison Is 
made between the profile loss coefficient data of Ainley and the data 
of Markov. The blade shapes from which the Markov data are derived 






z 




( 41 ) 
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are very similar to those of Fig. 25. It would appear that such blade 
shapes would have higher profile losses than those based on airfoil 
designs which usually are used for lower flow deflections. The data 
comparison shows this is not the case. 

Of special interest in the Markov data in Fig. 10 is the fact 
that an optimum solidity of (s/c) ^ .625 to .650 is evident. Higher 
solidities apparently block the flow thereby increasing the profile 
loss, and lower solidities encounter separation effects since the 
blades act more like individual airfoils. The low point in the curve 
corresponds to guiding the flow in a channel, retarding separation 
onset, but not blocking the flow. 

2. The tip clearance losses shown in Fig. 9 seem to give a 
large range of values depending on which curve is chosen. The sit- 
uation is perhaps more reasonable when the ratio of radial clearances 
to blade height is restricted to realistic values less than .03. 

The loss estimate from the Ainley equation is obviously outside the 
flow deflection range intended for that relation. The other curves 
for tip clearance loss for reaction blading could be based on dif- 
ferent degrees of reaction corresponding to the fan aspect of the 
curves, i.e., high reaction blading would use the highest curve. 

Some reasonable estimate could be made. The only relation for 
impulse blading available is that of Markov. The Meldahl relation 
is the only one which accounts for the trailing vortex inducing a 
loss even though the tip clearance is zero, but the dependence on 
blade aspect ratio appears unreasonably strong for this data as a 
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whole. For (h/b) * 1.00, the loss at (k/h) ■ .02 Is greater than the 
highest ordinate shown. 

3. The Markov data Is all that Is available for making an 
estimate of trailing edge thickness loss directly. The plots In 
Fig. 20 Indicate the trends which would be expected. 

The Alnley method of modifying the complete loss coefficient 
depending on the relative trailing edge thickness seems an unfor- 
tunate way to present the data. The system is workable but indirect. 
It would be more useful for comparison and design purposes to have 
the trailing edge thickness data presented directly. 

4. The situation regarding secondary flow loss coefficients is 
the most confused of all. The general equation by Vavra is not very 
precise but certainly Indicates the proper trend if the model pre- 
sented In Fig. 7 Is correct. None of the other presently available 
methods Investigated will even show the proper trends for high flow 
deflections. 

The Alnley relation plotted in Fig. 18 seems to have depended 
too heavily on stator data with an axial entering angle and low flow 
deflections in the rotor. Otherwise, the curves Indicate that for 
a given entering angle. Increasing the flow deflection decreases the 
loss. 

In Fig. 15 the relation derived from Meldahl and Betz Indicates 
the proper trend for low flow deflections only. Taking the plot for 
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• -45° as an example, it can be seen that for flow deflections 
greater than =• 90° the indication is that increasing the flow 

deflection decreases the loss. 

The Markov relation plotted in Fig. 11 suffers the same deficiency 
pointed out for the Ainley and Meldahl curves plus abnormally low 
secondary loss coefficients throughout for reasonable blade aspect 
ratios. The data is of no value unless used specifically with Markov 
relations for all losses. In Ref. 20, when sample loss estimations 
are carried out, the total loss coefficients compare well with the 
measured values. Apparently the profile, tip clearance and trailing 
edge loss coefficients with the blade height corrections used in the 
complete method counter the low secondary loss estimates. 

5. The relations now available can be used to make a quantitative 
comparison of the component losses. The turbine designed in the 
following section will be used as a model to show the comparisons. 

The physical dimensions of the stator and rotor are as listed in 
Table IV and these are used in the relations and plots so far developed. 
The numerical results are shown in Table V for the stator and Table VI 
for the rotor. 



- 33 - 



Table IV 

TURBINE PHYSICAL DIMENSIONS 



1 


Stator 


Rotor 


Entering angle 


0° 


69.60° 


Exiting angle 


75° 


-71.57° 


Blade height, h 


1.165 in 


1.465 in 


Trailing edge thickness, t^ 


.036 in 


.025 in 


Blade chord, c 


3.16 in 


1.355 in 


Blade width, b 


1.60 in 


1.375 in 


Spacing, s 


1.80 in 


.9415 in 


Blade thickness, t 


. 69 in 


.6485 in 


Tip clearance assumed, (k/h) 


o 

o 


.02 



Table V 

STATOR LOSS COEFFICIENT COMPARISON 



Author 


Profile 

Loss 


Tip 

Clearance 

Loss 


Trailing 

Edge 

Loss 


Secondary 

Loss 


Vavra, Ref.l 








.070 


Ainley, Ref. 23 


.067 








Ainley, Ref. 25 


.047 






.041 


Markov, Ref. 20 


.030 




.032 


<.005 


Meldahl, Ref. 21 




.106 




.022 


Stodola, Ref. 3 


.040 
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Table VI 

ROTOR LOSS COEFFICIENT COMPARISON 



Author 


Profile 

Loss 


Tip 

Clearance 

Loss 


Trailing 

Edge 

Loss 


Secondary 

Loss 


Vavra, Ref.l 








.231 


Ainley, Ref. 23 


~ .12 








Ainley, Ref. 23 


.148 


.13 




.168 


Markov , Re f . 20 


.068 


.025 


.032 


.025 






★ 






Meldahl, Ref. 21 




.223 




.2601 


Stodola, Ref. 3 




.05* 







* 

reaction 



As can be seen, the component losses presented in Tables V and VI 
can be combined selectively to obtain almost any desired answer. 

The loss investigation is disappointing in that no clear-cut 
answers are provided. The complete Ainley method, or the' complete 
Markov method, may produce a reasonable estimate but the correlation 
is poor at best for component losses. 

Any desire to find a way of basing a new design on component 
loss considerations is presently thwarted. A need is evident for 
an orderly testing program to find the magnitude and interaction 
of component losses. 
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IV. Design of Turbine 



As a further investigation into the characteristics of a turbine 
being used for maximum work output, a single stage turbine is designed 
for a slightly supersonic velocity leaving the stator and a relative 
Mach number of 0.8 as seen by the rotor. Then the overall pressure 
ratio is increased and the off-design performance evaluated. The 
number of calculations required are prohibitive for hand calculations 
so the equations necessary for the performance evaluation are also 
programmed in Fortran language for the CDC 1604 Digital Computer. 

A. Loss Coefficient Data 

A set of loss coefficients for the design were provided by 
Vavra from some unpublished data. These loss coefficients are shown 
in Figs. 21, 22 and 23 and will be used for the design and the per- 
formance evaluations. 

The abscissa of the plot of stator loss coefficients, is 

the Mach number for an isentroplc expansion from the total pressure 

at the inlet to the static pressure at the discharge of the blade 

row. This can be seen in the T-s diagram of the expansion through 

a turbine shown in Fig. 24 to be the isentroplc velocity V^, divided 

tn 

by the acoustic velocity based on T^^ When reduced, the result is 




( 42 ) 
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To properly determine the flow area through the stator, the 
loss coefficient in Fig. 20 noted "for flow rates" and shown in 
Fig. 24 as is used. This loss coefficient is intended to 

represent the losses encountered by the flow in the channel between 
blades, but is not intended to account for the mixing losses and 
separation losses encountered by the flow leaving the blade row. 
These additional losses are included in the higher loss coefficients 
which are applied to velocity determinations. These loss coef- 



Slmllar considerations are used for the loss coefficients through the 
rotor. ’ 

The loss coefficients are corrected for blade height according 
to the curves given in Fig. 23. 

B. Flow function 

For steady adiabatic flow the stagnation enthalpy remains con- 
stant along a streamline. For stators, the stagnation enthalpy is 



and in Fig. 24, to a different scale, would be represented by 
T^ and hj^ would be represented by Tj^. A similar situation pertains 
to a relative stagnation enthalpy defined as 



for the condition such that the peripheral velocity at the entrance 



flclents are noted "for velocities" in Fig. 21 and as 





(43) 




(44) 
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of the rotor equal to the peripheral velocity at the exit of the 
rotor. As shown in Appendix 1, the mass flow rate for a given set 
of inlet conditions can be expressed by 












yj) { pj 



( 45 ) 



where n is the poly tropic exponent, is the total pressure ahead 
of the stator blade row and p^ is the static pressure behind the 



blade row. 



As shown in Fig. 24, a relative total temperature and relative 
total pressure can be defined as 



Ze/ - 


t; ^ 






/ V 


^1 ~ 


^'ifr 



(46) 



(47) 



The mass flow rate through the rotor is then 



XiT' — 






(45a) 



Rearranging the above expressions for flow rate in non-dimensional 
terms defines a flow function. 






( 48 ) 
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and 



(p — /"A ('A ' j -" 

/g^x/z /»'-/Li/ij 1^/ 

It will be necessary to investigate whether this flow function can 
be satisfied for each set of conditions imposed across a blade row. 
The choked flow condition corresponds to the maximum value of P/p 
which the blade row can accomodate. The development of the critical 
values of the flow function and pressure ratio can also be found in 
Appendix I to be 




A computer program to do these computations is explained in section ! 
below, 

C. Geometry and Velocity Triangles 

The final design of the stator and rotor blade rows is as shown 
in Fig. 25 for the mean diameter. The velocity triangles of the 
average flow velocities at each station are as shown in Fig. 26. 

The sign convention chosen is that positive angles and positive 
peripheral velocity components point in the same direction as U, the 
rotor peripheral velocity. 



(48a) 



(49) 



(50) 
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It is convenient to define an **exit plane** and a **discharge 
plane** leaving each blade row as shown in Fig. 25 by the and @ 
designations. The design of the stator and rotor is such that the 
minimum cross-sectional area of the channel between blades occurs 
at the exit plane ^ . For maximum work output, the pressure drop 
across the stator is great enough to cause choked flow at the exit 
plane plus an additional expansion of the flow from the exit plane 
to the discharge plane. The conditions for such after-expansions 
are as shown in Appendix II based on the treatment in Ref, 26, 

In Ref. 1 a similar development is shown for a supersonic flow 
entering a blade row. Since the flow leaving the stator is super- 
sonic, it is possible under certain conditions for the rotor to see 
a supersonic flow also and proper account must be taken of the shock 
losses involved, 

D. Dimensionless Parameters 

As has already been mentioned, the number of variables to be 
handled is greatly reduced if the mass flow expression is non-dimen- 
sionalized to form a flow function depending only on pressure ratio, 
specific heat ratio and the polytropic exponent of the expansion 
process . 

In Ref, 27 can be found a complete development of dimensionless 
parameters and referred values for minimizing the number of variables 
necessarily handled in presenting the data for a turbine performance 
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analysis. Those used in the present analysis are referred rpm, over- 
all pressure ratio, referred flow rate, and a power coefficient. The 
necessary coefficient forms are computed from 



Referred RPM ■ RPM/ / T 

' o 


(51) 


Overall pressure ratio ■ 


(52) 


Referred flow rate ■ w V T /P 

1 o o 


(53) 


Power coefficient ■ HP/P "/"t 

of o 


(54) 



The referred flow rate and power coefficient are constant in this 

evaluation. It is useful to refer all the conditions at any point 

in the turbine to inlet conditions. For this reason, velocities are 

carried as V/ ^nd W/ and pressures as p/P and P„/P . 

9 0 9 o o Ro 

E. Design Computations 

In order to design the stator for choked flow at the exit plane, 
the flow function must be known for the imposed conditions. A com- 
puter program was written to compute the flow function for a given 
specific heat ratio, ^ , pressure ratio, P/p, and loss coefficient, 

As a by-product the program also computes the poly tropic loss coef- 
ficient, f ThU information is presented in Appendix V as Table III, 
Flow Function and Polytropic Loss Coefficients, for X from 1.25 to 
1.40, from 0.0 to 0.25 and P/p from 1.02 to 6.00. The particular 
flow function data to be applied here for ^ » 1.37 are also included 
as Fig. 27. The computer program is included in Appendix III. 
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The following data were assumed given for the design: 



Mass flow rate 



■ 100. Ibm/sec 



Total pressure at inlet * 1000 psia 
Total temperature at inlet * 1460°R 



Specific heat ratio 



- 1.37 



Gas constant, R 



= 421.5 ft Ibf/lbm °R 



Mean diameter of stator * 23.0 in 



Mean diameter of rotor 



- 23.0 in 



Design rpm 



- 13000 



The stator blade design was provided by Vavra from unpublished 
data and is as shown in Fig. 25. It is now necessary to select the 
blade height and number of blades so the choked flow condition at 
the exit is satisfied. 

A blade height of 1.16 Inches is taken as a first estimate. Then, 
from Fig. 21 the loss coefficient for flow areas at M^^ * 1.00 is 
found to be * .053. From Fig. 23 the blade height correction 

is IL »» 0.872 and so 



Using linear interpolation on Fig. 27 for this loss coefficient 




between the maximum values of 




found that 



^ » .65811 

^ cr 




cr 



- 1.86 
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Recalling the flow function equation (48), It Is found that 



A 0 

'throat 

A 

TttgOAT 




13.82 



13.82 

.65811 



21.0 in^ 



From Fig. 25 the measured minimum distance between blades at 
the exit plane is a » .448 in. and s « 1.8 in. is given. The number 
of blades is then 

^ ^ =* 40.2 

^ (55) 

In order to prevent harmonics in the wake pattern from the stator 
as seen by the rotor, a prime number of stator blades is chosen. 

In this case, z > 41 blades will suffice. For 41 blades, 

J. = ^ (23.0) - 1.763 

41 

so all measured values from Fig. 25 must be scaled by the scale 
factor. 

Scale Factor " * 0.9794 

1 . O 

The blade height must now be 

” * TwrfrTa qVq/\ * 1*165 in (56) 

SI za 41.0 (.448 x .9794) 

Recomputing with a new loss coefficient based on this new blade 

height makes no measurable difference In exit area required so the 
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scator row design Is now fixed. To simplify the iteration on ^ a 
computer program was written to find the choked flow values of 
(P/P)gr ^p* computer program is included in Appendix III 

and the computed data are presented as Table Vll in Appendix V. 



The flow conditions at the stator exit plane are determined 
from Fig. 24 as follows: 



Po 



1 - .538 

1.86 



^ = /^) *■ . .8458 

Po \ Poj 



^(IL -"{a) -.1470 

^ U rJ 
_^=7iiinir -m.52 

% I y-' % 

When a check is now made of the Mach number at the exit plane in 
order to check for choked flow it is found that 




- 0.97 



This computed Mach number is the average Mach number of the flow at 
the minimum area. The highest Mach number at this area is equal 
to 1.00 but it is averaged with the Mach number of Che flow in the 
boundary layer as accounted for by Che loss coefficient. In 
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Appendix IV the derivation of the Mach number for choked flow In a 
polytropic expansion Is presented. The "critical average" Mach 
number Is found to be 



1 /■h'^p(r-0 ] t'-/ < 

This condition will also be encountered later In the computer program 
for the flow through the rotor and proper account must be taken of 
it there. 

Using the relation from Ref. 20 for flow exit angle and the 
measured blade exit angle from Fig. 25 of “ 73.9° gives 

5= = .268 

= 75° 



For the rotor, the blade height Is permitted to be approximately 
0.30 Inches greater than the stator blade height In order to permit 
the streamlines to expand and offer no resistance to flow leaving the 
stator but not have dead areas above and below the useful flow area 
of the rotor blade. The rotor blade height will be taken then as 
h^2 “ 1-^65 in. 

In order to satisfy the requirement for the rotor to see a 
relative Mach number of 0.8, some after-expansion Is required. In 
order to evaluate the proper amount of after-expansion, the equations 
in Appendix II were coded for the computer to solve for any applied 
pressure ratio between the exit and discharge planes for a given ^ , 
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Mach number and blade angle. This computer progrsn is included in 
Appendix III and the solution for “IT • 1.37 and o(.^ • 75° over a 
range of pressure ratios is presented here as Fig. 28. The computer 
program can be used for any specific heat ratiot blade angle, Mach 

t t ^ 

number and pressure ratio range. 

After some iteration, it is determined that p. /P * 0.460 should 

l o 

give approximately the correct pressure ratio across the expansion 
area from the exit plane to the discharge plane. 



Pe Pe 



.460 

.538 



0.855 



Now from Fig. 28 it is found that /^/0( • 0.2° and ^ ^ .0008. The 
total loss coefficient to be applii4 is based on a blade height of 
1.165 Inches for the stator, fiiara is .1045 from Pig. 21 and 

■ .838 from Pig. 23 for velocity computations. Therefore 






Now Che velocity leaving the stator is determined with Fig. 24 
as reference 



_ 

z 


rtf H 


[lie)-”"---. 


To 


(i-¥) 


^ j m .1892 (.9117) - .1725 


ft ■] 


« 131.60 

^ r-/ ^ 
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To find the relative angle and relative velocity of the flow, It Is 
necessary to construct the velocity triangle as shown In Fig. 26. 



LL 

yr 






7T_ ^ 

72.0. 



«**i 



34.18 



75° - 0.2° • 74.8® 



Mmi _ V/ 74.8° ^ 

i^ft 

4^ 74.8° 

7t“ fn 

- U, 

n~ iT fT 

»v'< _ K., / 



34.5 



- 127.0 



■ 92.82 



60^ 



>— - -99.07 

/T^ frT 

As a check on the requirement of a relative Mach number of 0.8 as 
seen by the rotor 






99.07 

123.9 



0.80 



To design the rotor, several considerations must be kept In 
mind. The flow deflection should be as great as Is practicable for 
maximum work output. The trailing edge thickness will be arbitrarily 
selected as .025 Inches, and the spacing and blade design will be 
chosen so the minimum cross section area of the flow passage between 
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blades will occur at the "exit** plane. The rotor blade width is 
selected as 1.375 Inches and impulse conditions only will be used 
for unchoked flow, i.e. , p^^ « p^. 

Recalling that a relative total temperature and relative total 
pressure are defined as 



..9253 

7^ 7p ^ 7^ Ztf^ 




The flow function as applied to the rotor is 

(f> y ^ 

7m. fl / 

and this expression will be used to determine the exit area of the 
rotor. For the prescribed impulse conditions 




♦ 6955 
.4600 



1.515 



and the value to be used on the abscissa of the plots for rotor 
loss coefficients in Figs. 22 and 23 




34.18 

131.61 



■ .2598 



(46) 



(47) 



(48a) 
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With this value the rotor loss coefficient for flows is found to be 



• .819 (.192) • .1572 

Now using Fig. 27 for this pressure ratio and loss coefficient and 
using linear interpolation the flow function is 

^ • .58953 

and from eq. (A8a) above 

Aj • 19.1,3' 

Aj • 32.44 in^ 

An impulse rotor operating under conditions as these should have 
a ratio 

and so ’ 

s • 0.7 (1.375) • 0.962 

For the mean diameter of 23.0 Inches, the number of blades on the 
rotor is then 

*2 * • 75.1 blades 

Since an odd number of blades was chosen for the stator the rotor 
should have an even number of blades to reduce the possibility of 
resonance induced vibrations and the resulting metal fatigue. 
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Arbitrarily choosing 76 blades 




• .9415 In 



and the distance between blades, a, to obtain the required area of 
32.44 square Inches is 



The blade design itself is based on a circular arc contour of 
the concave face. The geometric relations of the design have been 
prograonned by Vavra fo^ computer solution. Using this computer 
program and requiring solutions for 74, 76, 78, 80 and 82 blades, 
the printout of the computer solution Is presented as Pig. 29. The 
values necessary to draw the blade design are shown in Pig. 30 and 
have been used to draw Pig. 25. 

The solution for 76 blades appears valid and so the other solutions 
are not used. The resulting |-|^| ■ 0.691. A check on this spacing can 
be made by Brillngs Rulei^ an early rule-of-thumb for steam turbines 
which is found in Ref. 3 to be 




- .273 In 



The exiting blade angle can be seen from Pig. 25 to be 
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( 57 ) 



It) • 2-5 *1” V * - 90 - 

|jj • 2.5 sin 36.86° • 1.5 

(f ) • 



Also Ref. 28 presents a relation for spacing and blade angles 
based on an aerodynamic load coefficient which is found to 

have a nearly constant value over a wide range of turbine designs. 
The relation is 





0.8 



(58) 



and so» using proper angles as defined in Ref. 28 for this case 

jil . ^ - .704 

(b( 2.0 18.43“ si„ 38.83° 



The relations by Briling and Zweifel show at least that the choice 
of 76 blades is reasonable for this design. The relation of Zweifel 
in particular is based on two-dimensional considerations and does not 
account for blade height. Some unpublished test data available to 
Vavra indicates the value of 0.691 is a good one for this size turbine. 

To complete the computations necessary to determine the power 
output of the turbine for a given set of conditions, the velocity 
triangle after the rotor is computed. 



- 51 - 



I o© 



T 



• 4tW' «• 




• •* • i' '»•■ 



li 



For no after-expansion, 



/2 



B2 



-71. 57 ”, and finding loss 



coefficients In the manner previously described, 



♦ 460 
.6955 



o \ .2701 , 



\ *" - .926 

^ ~ n\^J 

., 02.0 

z \zz 



.8275 



z 



Wz 

iT 



7e/ 

— ^ • .8543 



I r-/ 77 



• 84.39 



Using relations as before with Fig. 26 as reference. It Is found 



fT 



45.93 



and the power output, with the proper conversion factors is 



HP 



■ |j ''ul • "2''u2> 



(59) 



These computations are verified by the computer program. The first 
page of the printout of the computer program Is presented as Fig. 31. 
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V* Computer Program 



A. General Considerations 

The computer program for the evaluation of the turbine performance 
is Intended to be as general as possible within the framework of 
the turbine design method used. Any single stage turbine so designed 
can be evaluated by entering new physical dimensions, entrance con- 
ditions and loss coefficient data. The static pressures to be im- 
posed ahead of and behind the rotor must be programmed by the user 
within the main program. The loss data must be of the same format 
as Pigs. 21, 22 and 23. 

B. Main Program 

1. The main program is intended to make the major computations 
and comparisons and provide control. The subroutines are provided 
to do the lengthy iterative processes and repetitive calculations. 

The flow charts for the main program and the subroutines are in- 
cluded with the program listings found in Appendix III. 

The basic concept of the computer program is as follows: 

a. With the input data given, compute the design condition 
including the after-expansion at stator discharge. 

b. Reduce the static pressure through the impulse rotor in 
increments, computing velocity triangles and power at each pressure 
increment. The shock condition for supersonic flow entering the 
rotor is computed. Also the flow function through the rotor is 
checked for the choked condition. 
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c. When the choked condition through the rotor is attained, 
the static pressure at the rotor discharge is lowered in increments 
and the after-expansion condition computed. The resulting power and 
velocity and thermodynamic conditions are computed for each pressure 
increment. 

d. The process is terminated when the loss coefficient due 
to the rotor after-expansion is reduced to zero. 

The program will also terminate if the loss coefficient due to 
after-expansion out of the stator goes to zero or if the loss coef- 
ficient due to an entry shock into the rotor goes to zero. Any con- 
tinuation of the program beyond these limiting conditions would lead 
to conditions which violate the Second Law of Thermodynamics, as 
explained in Ref. 1. 

The loss coefficient reaches a maximum and then declines to zero 
or negative values as the Imposed pressure ratio across the control 
area is decreased in increments for the after-expansion case. The 
exact reason for this is not known. A check was made to see if the 
axial component of the Mach number of the flow leaving the blade 
row reached M ■ 1.00 when the loss coefficient reduced to zero. It 
is found this is not exactly the case but that the axial component 
of the Mach number is less than one. Further work is needed on this 
point. 



2. Some pertinent comments may be useful about the main decision 
points in the main program. The flow chart in Appendix III shows that 
the first check is for the magnitude of the static pressure between 
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the stator and rotor compared to the static pressure in the exit plane 
of the stator. For a stator designed as a converging-diverging nozzle, 
the Mach number at the exit plane may be greater than * 1.00. If 
the static pressure between the stator and rotor is greater than the 
static pressure at the stator exit plane, the shock condition is 
computed at the stator exit. If the pressures are equal, there is 
no flow deflection, and if the static pressure between the stator 
and rotor is lower than the exit plane pressure, the after-expansion 
case is computed. 

The next decision is based on the relative Mach number as seen 
by the rotor. For the supersonic case, the shock of the flow in the 
rotor entry is computed. For subsonic flow the component of relative 
velocity in the direction of the rotor blade entering angle is used 
to compute flow conditions. 

The flow function and critical pressure ratios are used to 
determine choked conditions in the rotor. As can be seen in Appendix IV 
the Mach number of the flow at the rotor exit plane is an average 
value, and the critical Mach number varies with the loss coefficient 
as shown there. Since the impulse conditions for the rotor are 
specified, some care is taken by an iteration process to find the 
Pj^ • P 2 condition for choked flow before the program continues. 

Once the choked condition in the rotor is attained, the con- 
ditions forward of the rotor exit plane remain fixed. The Mach number 
of the flow at the rotor exit plane is set at M„ "1.00 and the 
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static pressure behind the rotor is lowered in increments. The after- 
expansion behind the rotor is then computed until the limiting con- 
dition is reached. 

C. Subroutines 

1. The subroutine called TRNGL computes the velocity triangle 
based on rpm, absolute velocity and direction or based on rpm, 
relative velocity and direction. 

2. The subroutine called AFTER computes the deflection of the 
flow due to shock or expansion of the flow out of a blade row. The 
flow angle is incremented from the blade angle and the pressure ratio 
is computed and compared to the imposed value. The process con- 
tinues until the flow deflection is found which corresponds to the 
pressure ratio imposed. The detailed equations are derived in 
Appendix II, based on Ref. 26. 

3. The subroutine for computing the flow function for a 
pressure ratio imposed across the blade row is named CPHI. The 
polytropic loss coefficient corresponding to the design flow 

loss coefficient is determined by iteration. Then the flow function 
to satisfy the imposed conditions, critical flow function, and 
critical pressure ratio are computed from the equations developed 
in Appendix I. The pressure ratio was found to be more sensitive 
in fewer significant figures than the flow function itself so 
decisions in the main program are based on pressure ratio comparisons. 
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4. Subroutine BEFORE is used to compute rotor entry shocks based 
on the relations in Ref. 1 which in turn are based on considerations 
similar to those of Ref. 26 for after-expansion. Since the flow 
deflection is known, an iteration process is not required and the 
computations are straightforward. For cases which are not directly 
solvable by the basic equations, the component of the velocity in 

the direction of the blade angle is taken as the useful resultant. 

5. The CURVE subroutine is written as a means of finding a 
data point on a curve. The loss coefficient curves "for velocities" 
from Figs, 21, 22 and 23 are stored in the computer in the form of 
arrays. The subroutine writes the second-order equation of the 
appropriate three stored data points nearest to the input abscissa, 
and then finds the ordinate on the curve with the input abscissa. 

The subroutine is arranged in such a manner that the stored data 
must have equal abscissa increments and the first abscissa point 
must be zero. This means arbitrarily extending a curve such as 
Fig. 21 to M. ■ 0.0 whether the data at low values of M. will be 
used or not. 

D. Results of Computer Program 

1. The turbine design data given in the preceding section of 
this thesis were entered into the computer program and the resulting 
printout for the design condition only is presented as Fig. 31. The 
pressure ratio across the turbine was then increased in increments 
(by lowering p^^ and while holding constant) and the resulting 
computations are included as Table VIII in Appendix V. The rpm was 
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also decreased from 13000 to 9000 in 1000 rpm increments and the 
complete pressure range imposed at each rpm. These results are 
also included in Table VIII. 

2. The non-dimensional parameters previously discussed were 
then plotted in order to make cross plots which ultimately result 
in a portion of the turbine performance map. 

3. The data immediately available from the computer printout 
permits the direct plotting of Fig. 35, Turbine Efficiency vs. 
Pressure Ratio, with the rpm as a parameter. The anticipated gradual 
loss of efficiency with increasing pressure ratio is evident. Of 
interest is the change in slope of the curve at the point where 
rotor after-expansion commences. The indication is that the increase 
in velocity due to after-expansion causes the peripheral component 

of the velocity to increase, even though a deflection of the flow 
toward the axial direction also occurs. Therefore, the rate of loss 
of efficiency is decreased, although the change is very small. A 
reduction in rpm contributes to a large loss in efficiency at any 
given pressure ratio. 

4. In Figs. 36 and 37 the most important indication is as 

follows: Fig. 37 shows that for a fixed design at a pressure ratio 

of (^q/P 2 ) * 2.174, the increase in power for an increase in overall 
pressure ratio is direct and marked up to a pressure ratio of about 
4.5 or 5.0. Further large increases in pressure ratio net very 
small gains. In Fig. 36 the corollary is shown for efficiency. The 



- 58 - 



initial increase in pressure ratio and power output causes a certain 
loss of efficiency but the remaining increase in pressure ratio and 
power output causes a drastic loss of efficiency. 

An interesting comparison can now be made with the high head 
coefficient data from the first section of this thesis. In that 
section, each point on each curve represents a different turbine 
design. The conclusion was that if a required high power output 
transcends the need for high efficiency, then a very high head 
coefficient should be used for the turbine design. The initial 
increase of head coefficient from the optimum value causes a large 
loss of efficiency. Once this initial efficiency loss is accepted, 
further increase in head coefficient require little more loss of 
efficiency. Now in this section concerned with the off-design per- 
formance of a particular fixed turbine design, the conclusion is: 
additional power can be obtained at a higher pressure ratio than 
design but only in the initial increase of pressure ratio. The 
extremely high loss in efficiency is encountered if the last possible 
increment of power is demanded. 

5. The turbine performance map presented as Fig. 38 serves to 
indicate the pertinent trends for the parameters investigated. The 
lines of constant efficiency form only a portion of the elliptical 
curves they become over the complete range of referred rpm and 
power coefficient. The only conditions investigated in this thesis 
are at high power coefficients. The peak efficiency for this turbine 
would occur in the lower right hand corner of Fig. 38, i.e., in a 
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lower power coefficient range than shown here and in the referred 
rpm range of about 350. 



6. This turbine designed for high head coefficients can now be 
compared to the theoretical values predicted in the first section of 
this thesis. The head coefficient is defined in eq. (1), The 
Isentropic velocity is given by reference to the of Fig. 24 as 







(60) 



For pressures and temperatures at the design point 




Referring now to Table I at * 75^, r* * 0.0 and k, =» 16.0 the 

1 is 

values determined from the theoretical conditions can be compared 
to the turbine design. The comparison is shown in Table IX in the 
columns Theoretical #1 and Turbine Design. 



The ratio arbitrarily selected equal to one for the 

theoretical development. This was not a selected parameter for the 

turbine design and the discrepancy between the two values of V «/V , 

mz ml 

are the basis for the differences seen in Table IX. Since the ratio 

V «/V 1 was not required to be equal to one in the turbine design, 
mZ ml 

the flow is deflected more in the rotor. This results in higher 
absolute values for and Since the flow deflection is 

greater, the peripheral component of the absolute velocity (V^ 2 ^ 
greater and the power output and efficiency are increased. Otherwise, 
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the theoretical predictions are good. The from the 

turbine design can be entered into the equations and the corres- 
pondence should be more exact. This was done and the resulting data 
is the column headed Theoretical #2 in Table IX. 



Table IX 

COMPARISON OF THEORETICAL AND TURBINE DESIGN VALUES 





Theoretical #1* 

1 


Turbine Design 


r " 

Theoretical #2 




75° : 


74.8° (after-expansion) 


75° 


•k 




0 * 


0 


0 


•k 


R^/Rj 


1.0 * 

1 


1.0 


1 


k 


V/'^ml 


1.0 * 


.784 


.784 


k 




.95 * 


.955 


.95 


k 


^R 


.95 * 


1.00 


.95 


k 




.849 


.852 


.833 




'k 

r 


0.0 * 


0.0 


0.0 


k 


"'is 


16. * 


16.2 


16.2 


k 


7 


.60111 


.723 
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7. The curves of Fig. 26 are shown again in Fig. 39 with all 
values referred to the design point. All values are then shown as 

a percent of the design point values. Fig. 39 demonstrates the strong 
effect of rpm on power output and efficiency. At 707o design rpm, it 
is not possible to develop design power output. At about 907» 
design rpm, the design power can be developed at 957o of design 
efficiency. A maximum power of about 1257o of the design value is 
the maximum which can be developed at 907o design rpm. At the design 
rpm, 1257* design power can be developed with a loss of efficiency 
of 107«. When greater power increments are required, the efficiency 
loss becomes increasingly large. A 407o increase in power can be 
developed but at only 12 , 5 % of design efficiency. 

At design rpm, 1257o design power requires an«.overall pressure 
ratio of 3.5, and 1407o design power requires an overall pressure 
ratio of 8.5. 

8. As this thesis was being published, it was determined that 
the computer program finds the off-design performance of a turbine 
in which impulse conditions are maintained until choked. flow occurs 
in the rotor. A rotor blade height of 1.648 inches is required to 
satisfy such a condition. When this new rotor blade height is 
introduced into the computer program, the only direct influence is 
on the blade height correction factor for loss coefficients. The 
performance computations change is negligibly small. 
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VI. Conclusions and Recommendations 



A. Conclusions 

1. If a requirement for a high shaft power output turbine is 
encountered for which weight, size and space limitations preclude 
having several stages operating at peak efficiency, then a slightly 
higher head coefficient and slightly lower efficiency compromise 
design should not be made. Designing for power outputs higher than 
optimum should be based on extreme power output per stage. The 
initial loss of efficiency is great when a slightly higher than 
optimum head coefficient is considered. Not much more efficiency is 
lost by imposing extreme values of head coefficient. 

2. From the available literature, no straightforward method 

of computing losses is found which correlates well with the available 
data. This appears to be a fruitful area for future research and 
testing. 

a. The available data indicates the tip clearance loss is 
not zero even for impulse conditions and no clearance gap. The 
trailing vortices induced at the ends of the blades introduce 
velocity components which do no useful work. Clearance gaps within 
the annulus boundary layers should give tip clearance loss values 
about the same as zero clearance. As the clearance gap is increased, 
the loss increases. Also the loss increases as the degree of reaction 
is increased. A high degree of reaction means a high pressure drop 
across the blade row and such a pressure differential is conducive 
to high leakage across the blade tips. 
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The loss due to tip clearance is also affected by the relative 
size of the gap to the blade height. The effect of leakage across 
the tips of the blades on the secondary flow will make measurements 
difficult to separate into tip clearance and secondary flow losses. 

b. The profile loss is a function of blade thickness to 
chord ratio, spacing, flow deflection and probably general blade 
shape. This loss offers the most possibilities for simple and 
accurate measurements in a cascade test rig. 

c. The trailing edge thickness losses are probably in- 
sensitive below a certain ratio of trailing edge thickness and 
spacing, and then assume a real importance at higher values. These 
losses, too, should be relatively simple to measure. 

d. The secondary loss measurement offers the most difficulty. 
No known theory accurately predicts the loss. Indications are that 
the secondary loss is primarily dependent on the flow turning angle 
and ratio of blade height to the chord or axial width of the blade. 

The tip clearance effect is also of primary importance and will be 
difficult to divorce in the measurements. 

e. Once the four component losses mentioned above have a 
firm basis in measurements, it will be appropriate to consider 
Reynolds number and Mach number effects as they affect the losses in 
turbine blade rows. 

3. The computer program for the turbine performance analysis 
yields results as accurate as the loss data available. The program 
only requires inputs of loss coefficients in the specified form and 
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a design based on the flow function Cj) for area determinations. 

The analysis of the off-design performance of this particular 
turbine indicates that large gains in power output can be obtained 
at pressure ratios two or three times the design value if a 157o drop 
in efficiency can be accepted. Attempting to obtain the last 
possible increment of power, however, requires extreme values of 
pressure ratio and another large loss in efficiency. 

4. Further investigation into the after-expansion behind a 
blade row is needed. A cursory examination of the magnitude of the 
axial component of the Mach number of the flow leaving an after- 
expansion showed it ot be less than one when the loss coefficient 
decreased to zero. The exact reason for this loss coefficient be- 
havior is not known. 

B. Recommendations 

1. That an organized effort be made to measure and define 
component losses in blade rows in a form useful for design purposes. 

a. It would appear that trailing edge thickness effects 
and profile losses could be measured in cascade test rigs with a 
limited number of profiles, solidities and blade heights. Some 
correlation of these measurements could be made with Ainley and 
Markov. 

b. The tip clearance loss should be measured in a rotating 
turbine test stand. The direct effects would not be hard to measure 
if the effect on secondary loss was not present. Since this effect 
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is present, it will require study to obtain quantitative correlation 
and isolate the tip clearance loss from the secondary flow effects. 

c. The secondary loss will probably have to be deduced by 
subtracting the outer component losses from the overall efficiency 
as measured for the complete turbine blade row. No method is pres- 
ently known for measuring the secondary loss directly. 

2. Tests on the after-expansion behind a blade row should be 
made. Physical measurements may shed some light on the reason for 
the loss coefficient reaching a peak and then decreasing to negative 
values as the expansion is increased. Further theoretical study is 
also required, especially to consider three-dimensional effects. 

The testing program should include this variable. 
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•54052 


.52809 


.51576 


.53947 


•52696 


.51457 


• 50230 


• 52635 


.51389 


.50157 


.48937 


.51375 


.50135 


.48910 


•47700 


.50169 


.48936 


.47719 


•46518 


• 4901 6 


.47791 


•46582 


• 45390 


.4791 5 


.46698 


.45498 


.44317 


•46864 


.45655 


•44465 


.43293 


•45861 


•44660 


.43479 


•42318 


•44902 


.43710 


.42539 


•41388 


•43987 


•42803 


•41642 


• 40501 


•4311 1 


.41937 


.40784 


• 39654 


.42274 


.41108 


.39965 


.38845 


.41472 


•40315 


.39182 


•IlSIl 


.40704 


.39555 


.38431 


•37331 



• 200 



.225 



.250 



.17503 

.24228 

.29055 

.32865 

.36006 

.38663 

.40950 

.42942 

•44691 

.46238 

.47612 

.48838 

.49935 

.50920 

•51804 

.52600 

.53316 

.53962 

.54544 

.55068 

.55539 

.55963 

.56343 

.56684 

.56988 

.57259 

.57500 

.57712 

.57899 

.58061 

.58201 

.58321 

.58422 

.58505 

.58572 

.58624 

.58662 

.58687 

.58699 

.58700 

.58691 

.58673 

.58644 

.58608 

.58563 

.58511 

.58452 

.58387 

.58316 

.58239 

.57235 

.55965 

.54578 

.53153 

.51736 

.50351 

.49014 

.47730 

.46502 

.45331 

.44214 

.43151 

.42139 

.41175 

.40256 

.39379 




.17225 

.23840 

.28587 

.32330 

.35416 

.38025 

.40270 

.42223 

.43938 

.45453 

.46799 

.47998 

.49071 

.50033 

.50896 

.51673 

.52371 

.52999 

.53565 

. 54074 

.54532 

.54942 

.55310 

.55639 

. 55933 

.56193 

•56424 

.56627 

.56805 

.56959 

.57091 

.57204 

. 57297 

. 57374 

.57435 

.57481 

.57513 

.57532 

.57540 

.57537 

.57523 

.57500 

.57468 

.57427 

.57379 

. 57324 

.57262 

.57193 

.57119 

.57040 

•56016 

.54737 

.53347 

.51925 

.50513 

.49136 

.47808 

.46534 

.45317 

.44157 

.43052 

.42001 

.41001 

.40049 

.39141 

il94 




. 16942 
•23446 
.28111 
.31788 
.34817 
.37378 
.39579 
.41494 
.43174 
.44658 
.45975 
.47148 
.48196 
.49136 
.49978 
.50735 
.51416 
.52027 
.52577 
.53071 
.53515 
.53912 
•54268 
.54586 
.54868 
.55119 
.55340 
.55534 
.55703 
.55849 
.55974 
.56079 
.56166 
.56237 
.56291 
.56332 
.56359 
.56373 
.56376 
•56368 
.56350 
.56323 
. 56287 
.56242 
.56191 
.56133 
• 56067 
.55996 
.55919 
.55838 
.54796 
.53511 
.52120 
.50701 
.49297 
.47929 
•46611 
.45349 
•44144 
.42996 
•41904 
•40866 
•39878 
•38938 
.38044 
17191 

lip 

14863 
[41S6 
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HACH 


ALFA 


GAMMA 






1.000 


75 


r.370 


• 




PRAT 


PRAT1 


VELRAT1 


OALFA 


2A 


1.0000 


1.0000 


1.0C00 


-.0000 


.0000 


.9800 


.9800 


1.0146 


.0039 


.0000 


.9600 


.9600 


1.0292 


.0157 


.0000 


.9400 


.9401 


1.0438 


.0355 


.0001 


.9200 


.9201 


1.0583 


.0636 


.0001 


.9000 


.9001 


1.0729 


.1002 


.0003 


.8800 


.8799 


1.0876 


.1460 


.0005 


.8600 


.8601 


1.1021 


.2001 


.0008 


.8400 


.8401 


1.1167 


.2641 


.0011 


.6200 


.8201 


1.1313 


.3380 


.0016 


.8000 


.7999 


1.1461 


.4230 


.0021 


.7800 


.7600 


1.1606 


.5180 


.0028 


.7600 


.7599 


1.1753 


.6250 


.0036 


.7400 


.7399 


1.1900 


.7440 


.0045 


.7200 


.7200 


T.2045 


.8750 


.0055 


.7000 


.7001 


1.2191 


1.0200 . 


.0067 


.6800 


.6800 


1.2338 


1.1810 


.0080 


.6600 


.6600 


1.2485 


1.3570 


.0095 


.6400 


.6400 


1.2632 


1.5500 


.0112 


.6200 


.6199 


1.2780 


1.7620 


.0130. 


.6000 


.6000 


1.2928 


1.9930 


.0150 


.5800 


.5799 


1.3076 


2.2460 


.0171 


.5600 


.5599 


1.3225 


2.5230 


.0194 


.5400 


.5399 


1.3374 


2.8250 


.0219 


.5200 


.^5199 


1.3525 


3.1560 


.0246 


.5000 


.4999 


1.3676 


3.5180 


.0274 


T.IHEf 1 


MINUTES 


ANO 36 SECONDS 





F/G.23 

APTER-EXPANS/ON AT 137 AND 
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m- 



BLADING FOR ROTOR DESIGN 



DIAMETER 


S 


23.0000 


23.0000 


23.0000 


NUMBER OF BLADES 


; » 


. 74 


76 


78 


INLET ANGLE 


s 


20.4000 


20.4000 


20.4000 


EXIT ANGLE 


« 


18.4300 


18.4300 


18.4300 


AXIAL WIDTH 


• X 


1 .3750 


1.3750, 


1 .3750 


SPACING 


X. 


.9764 


.9507 


.9264 


EDGE THICKNESS 


X 


.0250. 


.0250 


.0250 


R 


X 


.7291 


.7291 


.7291 


XI 


X 


.2541 


.2541 


. .2541 


X2 


k 


.2305 


.2305 


.2305 


Y1 


X 


.6833 


.6833 


.6833 


Y2 


X 


.6917 


.6917 


.6917 


X3 


s 


.32 59 


. 3259 


.3259 


X4 


X 


.3096 


.3096 


.3096 


RS 


X 


.4276 


.4362 


.4443 


X5 


X 


.9821 


.9562 


.9317 


Y3 


X 


.7002 


.6998 


.6994 


Y4 


X 


.6748 


^.6752 


.6756 


X6 


X 


2.2681 


2.2681 . 


2.2681 


Y5 


X 


.7223 


.7223 


.7223 


Y6 


X 


.6527 


.6527 


,6527 


AREA 


X 


.661387 


.649473 


.637951 


DELTA 


X 


.9850 


.9850 


.9850 


CGtXS) 


X 


.9434 


.9361 


.9292 


CG(X7) 


X 


.9433 


.9360 


.9290 


CGI Y7) 


X 


.6996 


.6994 


.6993 


CG(Y8) 


X 


.6754 


.6756 


.6757 


I (X-X) 


* • 


06302621 . 


06249739 . 


06196594 


IIY-Y) 


* • 


03595965 . 


03399491 . 


03219183 


X(MAX) 1 


X 


.6378, 


.6305 


.6235 


YIMAX) 1 


X 


.6619 


.6619 


.6619 


X(MAX)2 


X 


.6259 


.•6186 


.61 16 


Y(MAX)2 


X 


.6876 


.6876 


.6876 


X(MAX)3 


X 


.4664 


.4564 


.4470 



TIMEt 0 MINUTES AND 33 SECONDS 



23.CC0C 

8C 

20.U000 
18.U300 
1 .375C 
.9032 
.0250 
.7291 
.2541 
.2305 
.6833 
.6917 
.3259 
-.3096 
. .4521 

.9084 
.6990 
.6760 
2.2681 
.7223 
.6527 
.626808 
.9850 
.9225 
.9223 
.6992 
.6758 
06143323 
03053410 
.6168 
.6619 
.6049 
.6876 
.4381 



23.0000 

82 

20.4000 

18.4300 

1.3750 

.8812 

.0250 

.7291 

.2541 

.2305 

.6833 

.6917 

.3259 

.3096 

.4594 

.8863 

.6986 

.6764 

2.2681 

.7223 

.6527 

.616030 

.9850 

.9160 

.9159 

.6991 

.6759 

06090044 

C2900724 

.6104 

.6619 

.5985 

.6876 

.4298 



FfG. Z9 
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TURBINE PERFORMANCE 



INPUT VALUES AND ^CONSTANTS COMPUTED 

FLOW RATE « 100*0 LBM/SEC PO » 1000* PSIA 

" >1 GAMMA « T.r“ 



REFERRED FLOW RATE 

Cl » 3.6195 

316*9026 
*000010 
136*3049 



3*1 



GAMMA 



,37 



TO s 
R « 



1460* 3 
421*50 



C4 



EXP1 « *2701 
EXP2 * 3*7027 
EXP3 « 2*7027 
RTO = 38*21 



RPM 



13000* 



PHYSICAL PARAMETERS 

REFERRED RPM 



MEAN OIA AT ROTOR INLET 
MEAN OIA AT ROTOR EXIT 



23*00 IN 
_ 23*00 IN 

STATOR DISCHARGE ANGLE 
ROTOR INLET ANGLE 
ROTOR DISCHARGE ANGLE 



STATOR BLADE HEIGHT 
ROTOR BLADE HEIGHT 
> 75*00 DEG 

* 69*60 DEG 

« -71*57 DEG 



LOSS COEFFICIENT THROUSH ROTOR (FLOWS! 
STATOR EXIT PLANE 



*1572 



340*226 

« 1*165 IN 
« 1*465 IN 



MIVEI 

1*00 


VE/RTO 

121*515 


PE/PO 

*5380 


TE/TO 
*85 30 




. 


STATOR 


DISCHARGE 


PLANE 








Mivn 

1*06 


Vl/RTO 

131*612 


Pl/PO 

*4600 


Tl/TO 
*82 75 


ALFA1 

74*78 


ZETA1 

*0883 


M(W1 1 
*80 


W1/RT0 

99*071 


PRl/PO 

*6955 


TRl/TO 
*92 53 


BETA1 

69*60 




IN ROTOR ENTRY 










M(W1) 

*80 


Wl/RTO 

99*071 


PRl/PO 

*6955 


TRI/TO 
*92 53 


BETAl 

69*60 




FLOW FUNCTION COMPUTATIONS 

PHI COMPUTED « *58953 
PRESSURE RATIO IMPOSED > 1*512 
ROTOR EXIT PLANE 


CRITICAL 


PHI CRITICAL 
PRESSURE RATIO 


M(WE) 

*73 


WE/RTO 

90*952 


PE2/P0 

*4600 


TE2/T0 

*8*29 






ROTOR 


DISCHARGE 


PLANE 








M(W21 

*67 


W2/RT0 

84*389 


P2/P0 

*4600 


T2/T0 

*6543 


BETA2 

-71*57 


ZETA2 

*2744 


M(V2) 

*42' 


V2/RT0 

53*105 


P2/P0 

*4600 


T2/T0 

*8343 


ALFA2 

-59*84 





*60837 

1*827 



RESULTANT POWER AND EFFICIENCY 

VU2(FT/SECI P2(P’SIA1 
-1754*49 460* 

OVERALL PRESSURE RATIO « 2*174 HP * 48682* 



VUK FT/SEC 1 
4852*61 



P2T(PSIA) 

519k 



T2 

1247* 



T2T 

1288* 



EFFICIENCY » *723 
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FIG. 34: 
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Appendix I. Development of Flow Function Relations 



In Fig. 32 is shown the adiabatic expansion process in a 

turbine from state 1 to state 2. 

The fluid proceeds from state 1 to state 2* and the isentropic 

enthalpy drop is represented by process The fluid proceeds 

from state 2* to state 2 and the isentropic enthalpy drop is process 

(2* -2 ). For the complete process the isentropic enthalpy drop is 

represented by process However, the sura of processes 

"k 

+ (2*-2 ) is greater than process (1-2^^) due to the so-called 
"reheat factor". The frictional effects raise the temperature, and 
therefore the enthalpy, of the fluid as it proceeds from state 1 to 
state 2*. This increase of enthalpy is then available to do work. 

As a result, it is necessary to make a distinction between the 
loss coefficient for the overall process of Fig. 14 and the loss 
coefficient f ^ for the polytropic process fro„ pressure (p) to 
(p-dp) on Fig. 32. 

The polytropic loss coefficient is 

i^p ■ 1 - /'p 



Looking at the expansion process in more detail in Fig. 32, 
the efficiency can also be seen to be 







( 1 . 2 ) 
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From this the relation along the polytropic line is 










= 7 






( / / 




2 !' 



and 






or 



^-r ^ -y-/ 

T >' 'IS' f 

Since, for an adiabatic process 

a constant, then / cC 

and for an isentropic process 

K' 

a constant, then / cC. 






Comparing these exponents, there is 
yr\ ^ / 






- 7 .(-^) ■ 



and 



/TL =r 



V 



'fpC'f-') ■+ I 



(1,3) 



(1.4) 



(1,5) 



( 1 , 6 ) 



(1,7) 
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V ^ 







To find the relation between the overall loss coefficient and the 
polytropic loss coefficient, eq. (1,5) can be written between two 
points on the polytropic expansion line as 



51 = 

T 



5 

r 



'h V-! 



and 



r-7-= Z 





% 




r-T 






= 7 ^ 




U) 




/ 





( 1 . 8 ) 



Solving equation (1,8) for and using 
the reduced result is 

f = 



-r-i 
-f \ 'r 









' - 

To define a flow function for the mass flow through the blade 
row, the expression for mass flow rate must be expressed in terms of 
the pressure ratio and the polytropic exponent. The mass flow is 



(1,9) 



^ 









( 1 , 10 ) 
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and following the procedure outlined in Appendix IV for finding the 
velocity corresponding to a pressure drop the result is 






Z 



T-/ 




(IV, 4) 



Along the poly tropic line, for an adiabatic process 

yrv 




* constant 



or 









/ 







( 1 , 11 ) 



When equations (IV, 4) and (1,11) are substituted into (1,10) the 
result can be arranged in non-dimensional form to define a flow 
function (j) to be 




jt 




( 1 , 12 ) 



This flow function reaches a maximum at the choked condition for the 
blade row. To find the corresponding pressure ratio for choked flow, 
the derivative of the flow function with respect to the pressure ratio 
will be set equal to zero and the result solved for the pressure ratio. 
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Then 












'Sr-/ 



-—/I 

/TV 



xrt ^ 



-h 



^'h ! 



yrr+J^ 

y*x 



-/ 






which reduces to 



_ ■23' 



Z-h'y^ 






^-1 






/n. /y^ 






When set equal to zero, the result is 



yt 



'tLZI 



--{—) 
^ \ ^-h / / 



and 






z 

/r>-hf 






(1,13) 



When the relation for the critical pressure ratio, eq. (1,13), 
is substituted into eq. (1,12) the critical flow function is obtained 
as 




(1,14) 
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< 
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Appendix II. After-expansion from a Blade Row 



To find relations to express the expansion of a sonic or super- 
sonic flow out of a blade row, the following model is employed. The 
cascade is assumed to have infinitely many blades of zero thickness 
arranged as in Fig. 33 below. 




Fig. 33 Relations for After-expansion from a Cascade 

A unit height perpendicular to the page is assumed. Also the uniform 
adiabatic flow is assumed to be frictionless at so that R is 
perpendicular to ® . The control area under consideration is 
bounded by the surfaces y and , i,e. , the "wall", the 

"exit" plane and the "discharge" plane. is sonic or supersonic. 
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I 



The continuity equation is used to find a relation between the 



densities , 






and so 



\ ^<y( 






(II, 1) 



The momentum equation. Ref. 1, ignoring shearing stresses is 

^ (II. 2) 



where 



and so 



= 






/A 






^ X. Ao( ^ '<U^ 



= 



/S 

~ X 



/y>^ — 



A 

— X 



A 



~^s 



Of particular interest are the relations in the i direction 

^ =■ ^ i- ^ ’C- ^ Cfi-:L o( 
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and using 



'^5 — rjz. 'X 






k: J- 



reduces to 

L 



/C^ Ao< 



7-/ - ^ 

A 



= O 



Since 



:t = ±- -w ^ = 

IT 



\/ 

/CL 






/ - /U)tt Ao< 

li 






/ 



(11,3) 



Since the process is assumed adiabatic, the total enthalpy remains 
constant through the control area, 

z z 

and 

-/ = c^r = ^jL 7- = _»LX- 
J~ r-/ ' 'T-/ ^ 

and substituting 

= jlA. ^ 

■T-/ A. r-/ fU 

' (11,4) 
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Rearranging and substituting in eq. (11,1) and using 




results in 




c( 



l-h 




(11,5) 



Then equating eqs. (11,3) and (11,5) and collecting terms results 
in a quadratic equation in velocity ratio. 




The loss coefficient associated with the flow through the 
control area can be defined to be 







based on the relations in Fig. 34. Using perfect gas relations and 
Fig. 34 as a reference, the following relations are obtained for the 
loss coefficient. 
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I 



# * > t 







I -^*T**^ 

w 
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i 





The Mach number of flow behind after-expansion is then, using Fig. 34 
as reference for relations, 



\Vj\ ^/(^/ 

- :^J— = 

y-i a.J' V-! ^ 




Further 



/ ^ rzi 



'-a.. 












^ / 4- IC-f 












'Vx 



Z 





/ V- i’-/ Mf' 

z ^ 




(11,9) 
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Substitute into (11,9) above 




This computation is included in the computer program but is not 
printed out. It is inserted for ready use if a check is needed. 



( 11 , 10 ) 
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Appendix III. 



MAIN PROGRAM FOR TURBINE PERFORMANCE ANALYSIS 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

SUBROUTINE AFTER 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

SUBROUTINE TRNGL 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

SUBROUTINE BEFORE 
FLOW CHART 
PROGRAM NAMES 
PROGRAM LISTING 

SUBROUTINE CURVE 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

SUBROUTINE CPHI 
FLOW CHART 
FORTRAN NAMES 
PROGRAM LISTING 

PROGRAMS TURPLOT AND TURBINE (HIGH HEAD COEFFICIENT DATA) 
FORTRAN NAMES 
PROGRAM TURPLOT LISTING 
PROGRAM TURBINE LISTING 

PROGRAM FANNO 3 (FLOW FUNCTION AND ZETA POLYTROPIC) 
FORTRAN NAMES 
PROGRAM LISTING 

PROGRAM FANNO 4 (CHOKED FLOW VALUES OF FANNO 3) 

FORTRAN NAMES 
PROGRAM LISTING 

PROGRAM FOR AFTER-EXPANSION COMPUTATIONS 
FORTRAN NAMES SAME AS "AFTER" ABOVE 
PROGRAM LISTING 
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li PROGRAM 
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PROGRAM PERFORM 





FORTRAN NAMES, 


EQUIVALENT SYMBOLS AND DEFINITIONS 


AE2 


^e2 


Area at exit plane of rotor 


AJ 


J 


Mechanical equivalent of heat 


ALFABl 




Undeflected flow angle of absolute velocity 
leaving stator 


ALFAl 




Flow angle of absolute velocity leaving 
stator 


ALFA2 


0(2 


Flow angle of absolute velocity leaving 
rotor 


AME 


M 

e 


Mach number of flow at exit plane of stator 


AMIS 


M. 

IS 


Mach number for isentropic expansion 


AMVl 


“vi 


Mach number of absolute flow leaving stator 


AMV2 


"V2 


Mach number of absolute flow leaving rotor 


AMWl 




Relative Mach number of flow entering rotor 


AMW2 


«W2 


Relative Mach number of flow leaving rotor 


AMWE 




Relative Mach number of flow at exit plane 
of rotor 


AMWEX 




Average Mach number of flow at exit plane 
of rotor for choked flow 


BETAl 




Flow angle of relative flow entering rotor 


BETA2 

BETABl 


j^2 

y^Bl 


Flow angle of relative flow leaving rotor 

Relative flow angle for zero incidence 
entering rotor 


BETAB2 


^B2 


Undeflected relative flow angle leaving 
rotor 


CP 


C 

n 


Specific heat at constant pressure 



P 
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Cl 




(Constant #1) = "/ R/g 


C2 




(Constant #2) - V(2 IT gR)/(1(' -1) 


C3 




(Constant #3) = (ir-l)/(2 T gR) 


C4 




(Constant #4) = V '2(‘gR 


DM1 


"ml 


Mean diameter of blade at rotor entrance 


DM2 


"m2 


Mean diameter of blade at rotor exit 


DRPM 


^RPM 


Increment of rpm 


DTE 


(AT^/T„) 


Temperature drop from inlet to exit plane 
of stator 


DTF 


^/T 

o 


Final temperature difference from total to 
static temperature after rotor 


DTR 




Temperature drop between relative total 
temperature ahead of rotor and static 
temperature ahead of rotor 


DTRAT 


(4t/T ) 
o 


Temperature drop from inlet to discharge 
plane of stator 


DTE 2 




Temperature drop from inlet to exit plane 
of rotor 


DT2 




Temperature drop from inlet to discharge 
plane of rotor 


EFF 


7 


Efficiency based on total pressure after 
ro tor 


EXPl 




Exponent yj^l = (1T“l)/'?r 


EXP2 




Exponent #2 = T/(T“1) 


EXP3 




Exponent #3 = 1/(T»1) 


G 


g 


Acceleration of gravity 


GAM 




Specific heat ratio 


HKR 




Blade height correction for rotor loss 
coefficient 


HKS 




Blade height correction for stator loss 
coefficient 
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« «fl 



HP 




Horsepower 


HR2 


^2 


Rotor blade height 


HSl 


"si 


Stator blade height 


NMAX 




Number of data points stored for a 
particular loss coefficient curve 


PE 




Static pressure at exit plane of stator 


PE 2 




Static pressure at exit plane of rotor 


PHIC 




Flow function computed with input "2T and ^ 


PHICR 




Critical flow function for choked flow in 
rotor 


PPCR 


(P/P)„ 


Critical pressure ratio for choked flow in 
ro tor 


PO 


P 

o 


Total pressure ahead of stator 


PPRAT 




Pressure ratio through rotor 


PRAT 


Pl/Pe 


Pressure ratio between stator exit plane 
and stator discharge plane 


PRATO 




Overall pressure ratio 


PRP 


\l^’’l 


Relative pressure ratio 


PRl 




Relative total pressure ahead of rotor 


PRIY 




Relative total pressure behind shock in 
rotor entry 


PI 


Pl^^o 


Static pressure between stator and rotor 


P2 


Pr/'o 


Static pressure after rotor 


P2T 


P2 


Static pressure after rotor 


R 




Gas constant 


RFR 


W V T /P 

/ o o 


Referred flow rate 


RPM 


RPM/ y T^ 


Referred rpm 
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tm* 
* I 

4v 



Ji\l 

^ • 



'is 



RPMA 



Actual value of rpm 



RTO 


fr 

' O 




TE 


T /T 
e o 


Temperature at exit plane of stator 


TEIS 


T /T 

eis o 


Isentropic temperature at exit plane of 
stator 


TE2 


T ,/T 
ez o 


Temperature at exit plane of rotor 


T1 




True interval between data points on loss 
coefficient curves 


TO 


T 

o 


Total temperature ahead of stator 


T1 


T./T 
1 o 


Static temperature between stator and rotor 


T2 


T^/I^ 


Static temperature after rotor 


T2T 


^2 


Static temperature after rotor 


TP2 


^2 


Total pressure after rotor 


TRAT 




Isentropic temperature ratio through stator 


TRATR 




Isentropic temperature ratio through rotor 


TRl 


\l'\ 


Relative total temperature ahead of rotor 


TT2 


^T2 


Total temperature after rotor 


TIY 


T * 
1 


Static temperature behind shock in rotor 
entry 


U1 


v/X 


Peripheral speed of blade at mean diameter 
at rotor entrance 


U2 


“ 2 ' ■/ ’'0 


Peripheral speed of blade at mean diameter 
at rotor discharge 


UV 


1— 1 
> 


Abscissa for rotor loss coefficient data 


VE 


V / '{t 
e I 0 


Absolute velocity of flow at exit plane of 
stator 


VM 


V / Vt 
m • 0 


Meridional component of absolute velocity 


VUl 


V - / lT T 
ul / 0 


Peripheral component of absolute velocity 
leaving stator 
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VU2 




Peripheral component of absolute velocity 
leaving rotor 


VI 




Absolute flow velocity leaving stator 


V2 


V-/ y T 

2 » o 


Absolute flow velocity leaving rotor 


WDOT 


w 


Mass flow rate 


WE 




Relative velocity at exit plane of rotor 


WM 


w / 

m ' o 


Meridional component of relative velocity 


QUl 


“ui/yx 


Peripheral component of relative velocity 
entering rotor 


WU2 


\ 2 ^{\ 


Peripheral component of relative velocity 
leaving rotor 


W1 


"i^ 


Relative flow velocity ahead of rotor 


WIY 


w */ y T 

1 ' o 


Relative flow velocity behind entry shock 
in rotor 


W2 


“2/ 


Relative flow velocity behind rotor 


YY 




Ordinate of loss coefficient curves 


ZET 




Loss coefficient for velocity through rotor 
including af ter*expansion 


ZETAA 


/a 


After-expansion loss coefficient 


ZETAB 




Loss coefficient due to rotor entry shock 


ZETAS 


fs 


Stator loss coefficient 


ZETAAl 


fal 


Loss coefficient for flows through stator 


ZETAA2 


'fal 


Loss coefficient for flows through rotor 


ZETAl 




Loss coefficient for velocities through 
stator 


ZETA2 


fa 


Loss coefficient for velocities through 
rotor without after-expansion 
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1 18 



PROGRAM PERFORM 
DIMENSION YY(4,53 
COMMON Cl ,C2,C3,C 
READ 281,P0,T0,GA 
READ 282,RPMA,DM1 
READ 283, AME,PE.,Z 
DO 118 K « 1,4 
READ 284, J, (NMAX( 
NMAX = NMAX(J) 
READ 285 (YY(J,I) 
CONTINUE 
AJ = 778. 17 
G = 32.174 
RTO = SQRTF(TO) 
RFR = ■ WDOT»RTO/PD 
Cl=SQRTF(R/32. 174 
C2=SQRTF(2.*GAM*3 
C3=(GAM-l.)/(2.*3 
C4=SQRTF(GAM*32.1 
EXP1 = (GAM-1 . )/GA'1 
EXP2= (GAM)/ (GAM-1 
EXP3=1. /(GAM-1. ) 
CP = R*GAM/(AJ*(3 
PE = PE/PO 
TEIS = PE**EXPV 
DTE = ( l.-TEIS)*( 
TE - l.-DTE 
VE = C2*SQRTF(DTE 
■CALL CURVE ( 1 ,HS1 
CALL CURVE (2,HR2 
DO 200 J = 1,5 
DRPM = 1000.*FL0A 
RPMA = 13000. -OR? 
PRINT 250 

251 , WOOT,PD 

252, RFR, GAN 

253, C1,EXP1 

254, C2, EXP2 
25b,C3,EXP3 
256, C4, RTO 
RPMA/RTO 
257 

258, RPMA, R? 

259, DM1 ,HS1 

260, DM2, HR? 

261 , ALFABl 

262, BETAB1 

263, BETAB2 

264, ZETAA2 
= ALFAB1*3 



) ,NMAX(4 ),TI (4) 

4,GAM,EXP1 ,EXP2,EXP3,P0,T0,R 
M,R.WOOT . 

,DM2.HS1 ,HR2,ALFAB1,BETAB1 ,BETAB2 ,AE2 
ETAAl ,ZETAA2,AMIS 

J),TI(J)) 

,1=1, NMAX) 



1 1 1 



100 



101 



103 



PRINT- 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
RPM = 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
ALFABl 
BETABI = BETAB1*3 
BETAB2 = 6ETA62*3 
PRINT 265 
PRINT 266,AME,VE, 
L = 2 

DO 193 I = 1,30 
PI = 470. 

PI = Pl-FLOATF(I) 
PI = Pl/PO 
CONTINUE 
PRAT=P1/PE 
P2 = PI 

IF (Pl-PE) 100,13 
CALL AFTER (AME,P 
IF (ZETAA) 116,11 
GO TO 105 
CONTINUE 
ALFAl =ALFAB1 
ZETAA=0.0 
CONTI NUE 

CALL CURVE (3, AMI 
ZETAl = ZETAS*HKS 
TRAT = P1»*EXP1 



) 

2. 174*R/ (GAM-1 , ) ) 
AM*32. 174»R) 

74*R) 

. ) 

AM-1 . ) ) 

l.-ZETAAl) 

) 

,HKS,TI ,NMAX,YY) 
,HKR,TI ,NMAX, YY) 

TF(J-l) 

M 

tTO 

,R 



M 



.1415927/180. 
.1415927/180, 
. 1415927/180. 

PE,TE 



* 10 . 



1,100 

RAT, ALFABl, ALFAl , ZETAA) 
6, 103 



S,ZETAS,TI,NMAX,YY) 

+ZETAA 
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OTRAT = ( 1 ,-TRAT) 
T1 = l.-DTRAT 
VI = C2*SQRTF(0T^ 
AMVl ■ V1/(CU*S0^ 
IF (L-3) 119,120, 

120 PRINT 290 
L = 1 

119 CONTINUE 
PRINT 267 
AL = ALFAl*180./3 
PRINT 291 , AMVl, VI 
CALL TRNGL ( 1 ,V1, 
0TR=C3*W1*W1 
TR1=T1+0TR 
PRP = (TRl/Tl )**E 
PRl = P1»PRP 
AMWl = Wl/(C4*SQR 
PRINT 269 
BE = BETAl*180./3 
PRINT 260, AMWl, W1 
IF(AMW1-1.0)10U,1 

104 W1=W1 »C0SF(BETAB1 
PRINT 271 
0TR=C3*W1*W1 
TR1=T 1+OTR 

PRP = (TRl/Tl )**E 
PRl = P1*PRP 
AMWl = W1/(C4«S0R 
BETAl = BETABl 
BE = BETAl«180./3 
PRINT 268, AMWl, W1 
GO TO 106 

105 CONTINUE 

CALL BEFORE (AMWl 
PRINT 270 
W1 = WIY 

AMWl = W1/(C4*SQR 
PRl = PRIY 
BETAl = BETABl 
BE = BETAl*180./3 
PRINT 268, AMWl, W1 
IF (ZETAB) 114,13 

114 PRINT 288 
GO TO 198 

106 CONTINUE 
UV = Ul/Vl 

CALL CURVE 

ZETA2 = ZETAR*HKR 
PPRAT=PR1/P2 
TRATR = l./(PPRAT 
DTE2 = (TRl-TRATR 

, TE2 = TR1-DTE2 
WE = C2*SQRTF(DTE 
AMWE = WE/(C4»SQ^ 
PRINT 272 
CALL CPHI (ZETAA? 
PRINT 273,PHIC,PH 
PRINT 274,PPRAT,P 
IF (PPRAT-PPCR) 1 

107 CONTINUE 

' DELTA = ABSF(PPRA 
IF (DELTA-. 0005)1 

122 PE2 = P2 
PRINT 2-75 
PRINT 266, AMWE, W= 
PRINT 276 
DT2=( TR1-TRATR*TR 
T2 = TR1-0T2 
W2=C2*SQRTF(DT2) 
AMW2 = W2/(C4*S0R 
BETA2=BETAB2 • 
BE2= BETA2*180./3 
PRINT 291,AMW2,W2 



*( l.-ZETAl ) 
AT) 

TF(T1 ) ) 

120 



.1415927 

,P1,T1,AL,ZETA1 

DM1 , ALFA 1,RPM,VM,VU1 , WUl , U1 , Wl , BE T A 1 ) 



XP2 

TF(T1 ) ) 

.1415927 
,PR1,TR1 ,BE 
04,105 
-BETAl) 



XP2 

TF(TD) 

.141 5927 
,PR1,TR1,8E 

, BETABl ,BETA1 ,P1 ,W1,TR1, Tl, WIY, PRIY, Til 
TF(TIY)) 



.1415927 
,PR1 ,TR1 ,BE 
6,106 



(4,UV, ZETAR,TI,NMAX,YY) 



**EXP1 ) 

•TRl ) *( 1 .-ZETAA2) 

2 ) 

TF(TE2) ) 

,PPRAT,PHIC,PHICR,PPCR) 

ICR 

PC R 

07, 1 15,108 

T-PPCR) 

15,115, 122 



,PE2,TE2 
1 )*( 1.-ZETA2) 

TF(T2) ) 

. 1415927 

,P2,T2,BE2,ZETA2 



,ZETAB) 



. 126 - 



108 
1 10 

121 

115 



299 



112 

123 

109 



1 16 
1 17 



113 



CALL TRNGL ( 2 , W2, DM2 , BET A2 , RPM, WM , WU2 t VU2 , U2, V2 , AL FA2 ) 
PRINT 277 

AMV2 = V2/(C4*SQRTF(T2M 
AL2= ALFA2*190./3.1U15927 
PRINT 268,AMV2,V2,P2,T2, AL2, 

GO TO 113 
CONTINUE 
PRINT 286 
L = 6 

P2 = P2 + .00'01 
PPRAT = PR1/P2 

CALL CPHI (ZETAA2, PPRAT, PHIC,PHICR,PPCR) 

DELTA = ABSF( PPRAT-PPCR) 

IF (DELTA-. 0005)1 21 ,121, 1 10 
PI = P2 . 

GO TO 111 
CONTINUE 
AMWE •= 1.00 
PE2 = P2 

PPRAT = PR1/PE2 . 

TRATR = 1 ./(PPRAT**EXP1 ) 

0TE2 = ( TR1-TRATR*TR1 )*( 1 .-ZETAA2 ) 

TE2 = TR1-0TE2 
WE = C2*SQRTF(DTE2) 

AMWEX= WE/(C4*SQRTF(TE2) ) 

FORMAT (100X,F10.3) 

PRINT 275 

PRINT 266, AMWE,WE |PE2,TE2 
DT2={ TR1-TRATR*TR1 )*( 1.-ZETA2) 

T2 = TR1-DT2 
W2=C2*SQRTF(DT2) 

BETA2=BETAB2 

CALL TRNGL ( 2 , W2, DM2, BETA2, RPM, WM, WU2 , VU2 , U2 , V2 , AL FA2 ) 
AMW2 = W2/( C4*SQRTF(T2) ) 

AMV2 = V2/(C4*S0RTF(T2) ) 

PRINT 276 

6E2= BETA2*180./3. 1415927 
PRINT 291,AMW2,W2,P2,T2,BE2,ZETA2 
PRINT 277 

AL2= ALFA2*180./3. 1415927 
PRINT 268,AMV2,V2 ,P2,T2,AL2 
L = 3 
GO TO 113 
P2=P2-.01 

IF (L-5) 109, 123, 123- 

PRINT 290 

L = 1 

CONTINUE 

PRINT 287 

PRAT=P2/PE2 

CALL AFTER ( AMWE, PRAT ,BETAB2 , BETA2 , ZETAA ) 

IF (ZETAA) 116, lU, 117 

PRINT 289 

GO TO 199 

CONTINUE 

PPRAT = PR1/P2 

TRATR = 1 ./(PPRAT**EXP1 ) 

DT2=( TR1-TRATR*TR 1 )*( 1 . - ( ZETA2+ZETAA ) ) 

T2 = TR1-DT2 
W2=C2»SQRTF(DT2) 

CALL TRNGL ( 2 , W2, DM2 , BET A2 , RP M, WM, WU2 , VU2 , U2, V2 , AL FA2 ) 
AMW2 = W2/(C4*S0RTF(T2) ) 

AMV2 = V2/(C4*SQRTF(T2)) 

PRINT 276 

BE2= BETA2* 180. /3 . 1415927 

ZET = ZETA2+ZETAA 

PRINT 291 ,AMW2,W2 ,P2,T2, BE2,ZET 

PRINT 277 

AL2= ALFA2*180./3 . 1415927 
PRINT 268,AMV2,V2 ,P2,T2, AL2 
CONTI NUE 
P2T = P^*PO 
T2T = T2»T0 



DTF = C3*V2»V2 
TT2 = T2+0TF 

TP2 = P2T*( (TT2/r 2)**EXP2) 
PRATO * 1./P2 
TT2 * TT2»T0 
OHI = TO-TT2 

DH2 = T0*(1.-((TP2/P0)**EXPD) 
EFF = DH1/DH2 
VU1T = VUl*RTO 
VU2T * VU2*RT0 



HP = RTO*WOOT«l.+ 14»(VUlT.*ur-VU2T*U2)/(G*AJ) 

PRINT 273 

PRINT 279, VU1 T, VJ2T , P2T , TP2, T2T , TT2 
PRINT 280, PRATO, HP, EFF 
L = L + 1 

IF (AHWE-1.00) 198,112,112 

2500F0RMAT ( 1 Hi ,2 IX 19HTURBI NE PERFORMANCE// 1 3X35H INPUT 
ITANTS COMPUTED /) 

2510F0RMATI 13H FLOW RATE = F6.1, 17H LBM/SEC PO = F6.0, 
ITO = F6.0, 2H R ) 

2520F0RMAT (22H REFERRED FLOW RATE = F5. 3 ,4X8HGAMMA = 

/) 

= F9.4,21 X7HEXP1 
= F9.4,21X7HEXP2 
= F9.6,21X7HEXP3 
= F9.4,21X7HRT0 
PHYSICAL PARAMETERS 



VALUES AND CONS 



13H PSIA. 



1 F6. 2 

253 FORMAT 
FORMAT 
FORMAT 
FORMAT 



F4.2,7X5HR 



254 

255 

256 



( 7X5HC1 
(7X5HC2 
( 7X5HC3 
( 7X5HC4 



258 



FORMAT! /20X20H 
FORMAT! 7H.RPM 



F6.4) 
F6.4) 
F6.4) 
F6.2) 
/) 



3/) 

IN 



= F5.2,27H IN ROTOR BLADE HEIG 



262 

263 

264 

265 

266 

2670F0RMAT 
1T1/T0 

268 FORMAT 

269 FORMAT 
2700FORMATI23H 

10 TRl/TO 
2710FORMAT! 15H 
lETAl ) 

272 F0RMATI27H 



) 

) 

./) 



=7.0,23X15HREFERRE0 RPM = F12 
2590F0RMATI27H MEAN OIA AT ROTOR INLET = F5.2,27H 
1GHT = F5.3,3H IN ) 

2600F0RMATI27H MEAN DIA AT ROTOR EXIT 
IHT = F5.3,3H IN ) 

261 FORMAT! 12X26H STATOR DISCHARGE ANGLE = 

FORMAT! 12X26H ROI OR INLET ANGLE 
FORMAT! 12X26H ROIOR DISCHARGE ANGLE = 

FORMAT ! 7X41HL0SS COEFFICIENT THROUGH 
F0RMATM8H STATOR EXIT PLANE//32H M!VE) 

FORMAT !F6.2,F10. 3,2F8.4/1 

!23H STATOR DISCHARGE PLANE//49H M!V1) 

ALFAl ZETAl ) 

!F6.2,F10. 3,2F8.4,F8.2/) 

! 41H M!W1 ) W1/RT0 PRl/PO TRl/TO BETAl 
IN ROrOR BEHIND ENTRY SH0CK//41H M!W1) 

BETAl ) 

IN ROrOR ENTRY//41H M!Wl) Wl/RTO PRl/PO 



STATOR BLADE HEI 



F6.2,4H DES 
F6.2,4H DEG 
F6.2,4H DEG 
ROTOR ! FLOWS) 
VE/RTO 



= F7.4//) 
PE/PO TE/TO) 



/ 1/RTO 



P1/P0 



) 



Wl/RTO PRl/P 



TRl/TO 



B 



FLOW FUNCTION COMPUTATIONS) 

273 FORMAT I10X16H PHI COMPUTED = F6. 5 , 1 3X1 5HPH I .CRITICAL = F6.5) 
2740F0RMAT! 26H PRESSURE RATIO IMPOSED = F6.3,28H CRITICAL PRESSURE RA 
ITIO = F6.3) 

275 FORMAT ! 1 7H ROTOR EXIT PLANE//32H M!WE) WE/RTO PE2/P0 TE2/T0) 
2760F0PMAT !22H ROTOR DISCHARGE PLANE//49H M!W2) W2/RT0 P2/P0 T 
12/TO BETA2 ZETA2 ) ' 

277 F0RMATI41H M!V2) V2/RT0 P2/P0 T2/T0’ ALFA2 ) 

2780F0RMAT ! /15X31H RESULTANT POWER AND EFFICIENCY // 

1 3X1 1HVU1 ! FT/SE: ) ,2X1 1HVU2! FT/SEC) ,3X8HP2 ! PS I A) , 2 X9HP2T ! PSI A) , 3X 

2 2HT2,5X3HT2T) 



279 FORMAT 

280 FORMAT 
ICIENCY 

281 FORMAT 

282 FORMAT 

283 FORMAT 

284 FORMAT! 

285 FORMAT 

286 F0RMAT!/48H 

287 FORMAT !21H 

288 FORMAT 

289 FORMAT 

290 FORMAT 

291 FORMAT 

198 CONTINUE 

199 CONTINUE 



!F12.2,F15 .2,F1 1 .0,2F10.0,F7.0/) 

! 5X25HCVERALL PRESSURE RATIO = F6.3,2X5HHP = F6.0,15H 
= F4.3///) • 

I5F10.0) 

!9F8.0) 

I5F10.0) 

I3,I6,F10. 0) ■ 

!8F9.0)_ 

CANNOT SATISFY FLOW FUNCTION, INCREASE PI /) 
AFTEREXPANSION) 

IS NEGATIVE. ) 

IS NEGATIVE, ) 



EFFI 



ROTOR 

ROTOR 

ZETAB 

ZETAA 



! 1 9H 
! 1 9H 
! IHl ) 

!F6.2,F10. 3,2F8.4,F8.2,F8.4/) 
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200 



ALFABl = 
DETAB1 * 
BETAB2 = 
CONTINUE 
END 



ALFABl*! 80./3. 1U15927 
BETAB1*1 80./3.1U15927 
BETAB2*! 80 . /3. 1 4 1 592 7 
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SUBROUTINE AFTER 





FORTRAN NAMES 


, EQUIVALENT SYMBOLS AND DEFINITIONS 


ALFA 

COSRAT 




Undeflected flow ang(gfl' (blade angle) 
Cosine ratio of blade angle to flow angle 


DALFA 


Ao< 


c/g “ (degrees) 


DALFAR 


AoC 


</b " (radius) 


PRAT 


<Pd^Pe> 


Imposed pressure ratio 


PRATI 


<Pd/Pe>l 


Computed pressure ratio 


VD 


''d 


Velocity at discharge plane 


VELRAT 


''d'\ 


Velocity ratio of flow at discharge plane 
to exit plane 


ZA 




Loss coefficient due to after-expansion 
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SUBROUTINE AFTER ( AME , PRAT , AL FAS, ALFA 1 , ZA ) 

COMMON Cl ,C2,C3,:U,GAM,EXP1,EXP2,EXP3,P0,T0,R 
ALFA = A8SF ( ALFAS ) 

£E= ( 2./( (GAM-1.) *AME*AME) ) 

DALFA = 0.0 
M « 0 

400 M = M+1 

401 OALFA = (-. 1 ) »*(M-1 l+DALFA 
DALFAR = DALFA*3. 1415927/180. 

COSRAT = COSF { ALF A) /COSF{ ALFA-DALFAR) 

AA = GAM*AME*AME*C0SF(DALFAR)-C0SRAT*(GAM-1 . )*AME»AME/2. 

BB = -GAM*AME*AM=-1 . 

CC = COSRAT*( l.+( (GAM-1. )*AME*AME/2.) ) 

DO = BB*BB-4.*AA»CC 
IF (DO) ’401,402 , 402 

402 VELRAT = ( -BB+SQ^ TF ( DO) ) / ( 2 . * AA ) 

OPR ATI = (COSRAT*( 1 . + ( (GAM-1 . ) * AME* AME/2. ) * ( 1 . -VELR AT*VELR AT ) ) ) 
1 /VELRAT 

IF (PRATl) 410,411,411 

410 PRINT 409, N, AME, PRAT, ALFA, DALFAft 
409 FORMAT (72X, I3,4F 10.5 ) 

411 CONTINUE . * 

ZA = l.-(VELRAT*VELRAT)'/(l.+EE-EE*(PRATl**EXPin 

VD = VE*VELRAT 

DIFF = ABSF(PRATl-PRAT) 

IF (DIFF-. 00001 ) 408,408,403 

403 IF (PRATl-PRAT) 404,408,401 

404 M=M+1 

405 DALFA = ( - . 1 ) »* ( 'I -1 ) +OALF A 
DALFAR = DALFA*3. 1415927/180. 

COSRAT = C0SF(ALFA)/C0SF( ALFA-DALFAR) 

AA = GAM*AME*AME*COSF (DALFAR)-C0SRAT*(GAM-1. )*AME*AME/2. 

BB = -GAM*AME*AME-1 . 

CC = COSRAT*( l.+( (GAM-1. )*AME*AME/2.) ) 

DD = BB*BB-4.*AA*CC 
IF(DD) 400,406,436 

406 VELRAT = ( -BB +SQ^ TF ( DD ) ) / ( 2 . * AA ) 

OPRATl = (COSRAT*( l.+ ( (GAM-1 . )*AME*AME/2. )*(1.-VELRAT*VELRAT) ) ) 
1 /VELRAT 

IF (PRATl) 412,413,413 

412 PRINT 409, N, AME, PRAT, ALFA, DALFAR 

413 CONTINUE 

ZA = l.-(VELRAT*i/ELRAT)/(l. + EE-EE*(PRATl**EXPD) 

VD = VE»VELRAT 

DIFF = ABSF(PRATl-PRAT) 

IF (DIFF-. 00001 ) 408,408,407 

407 IF(PRAT-PRAT1 ) 430,408,405 

408 CONTINUE 

IF (ALFAS) 414,415,415 

414 ALFAl = ALFAS+OALFAR 
RETURN 

415 ALFAl = ALFAS - OALFAR 
END 
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SUBROUTINE TRNGL 



FORTRAN NAMES. EQUIVALENT SYMBOLS AND DEFINITIONS 



A 

B 

N 

U 

VM 

VUl 

WUl 

W1 



m 



ul 



W 



ul 



W, 



Input angle of flow 

Output flow angle (If A = , B 

B = ^ 




and if A 




Value selects plus or minus SIGN 
Peripheral velocity of blade at mean diameter 
Meridional component of absolute velocity 
Peripheral component of absolute velocity 
Peripheral component of relative velocity 
Relative velocity 
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SUBROUTINE TRNGL ( N , V 1 , DM , A, RPM, VM, VU 1 t WUl , Ut W1 » fl) 
COMMON Cl .C2,C3,CU,GAM,EXP1,EXP2,EXP3,P0,T0,R 
SIGN = (-1)»*N 
U = 3 . 1415927 *RP' 1 *DM/r 20 . 

VM = VI »COSF( A) ■ 

IF (A) 701,702,732 

701 VUl = -V1*SINF(A) 

WUl = VU1+SIGN»U ■ . 

B = ATA71F ( WUl/VM) 

W1 = WUl/SINFIB) 

702 CONTINUE * • 

VUl = V1-*SINF(A) 

WUl = VU1+SIGN*U 
B = ATANFIWUl/VM) 

W1 = WUl/SINFIB) 

END 
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SUBROUTINE BEFORE 





FORTRAN NAMES 


, EQUIVALENT SYMBOLS AND DEFINITIONS 


A 




Interim computation for velocity ratio 


ALFA 


"^R 


Blade angle entering rotor 


ALFAI 




Flow angle ahead of rotor 


AMD 




Mach number of flow after shock 


B 




Interim computation for velocity ratio 


cc 




Cosine ratio of angles involved 


D 




2 

A * B, used to determine whether shock should 
be computed, or the component of the velocity 
in direction of blade be used. 


DALFA 




= o('^ - 


DTR 


(T^l-ip 


Temperature difference between relative total 
temperature and static temperature 


FF 




Interim computation to find Mach number after 
shock 


GG 




Interim computation to find Mach number after 
shock 


PP 


Pj/Pl* 


Static pressure ratio through shock 


PRP 


‘■ri'P! 


Pressure ratio through shock 


PRIY 




Relative total pressure behind shock 


PlY 




Static pressure behind shock 


W 




Velocity ratio of flow ahead of and behind 
shock 


ZETAB 


' B 


Loss coefficient of flow due to entry shock 
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SUBROUTINE BEFORE ( AM 1 , ALFA , ALF A1 , P 1 , W 1 , TR 1 , T 1 , W 1 Y , PRl Y , T1 Y, ZET AB ) 

COMMON C 1 ,C2,C3,C4, GAM, EXP 1 , EXP2, EXP3,PO,TO.,R 

DALFA = ALFAI-ALFA 

c: = COSF (ALFAl )/ C0SF( ALFA) 

A * ( 1. /{ GAM+1 . n *{( 1 ./(CC»AM1 »AMn )+GAM*COSF (DALFA) ) 

B = -n .-(2./(GAM+l.) )*( ( l./( AM1*AM1) )+GAM) ) 

0 = A*A-R 

IF (D) 602,601,63] 

602 WlY = Wl-»COSF( ALFAI-ALFA) 

2ETA3 = 0.0 

DTR = C3*W1Y*W1Y 

T1 Y = TRl-DTR 

PRP = (TR1/T1Y)»*EXP2 

PRIY = P1Y*PRP 

RETURN 

601 CONTINUE 

VV = A-SORTF( A*A-B) 

PP=CC*( ( 1 . + (GAM-1 . ) *AMl*AMl/2. )*( 1 ./VV)-( (GAM-1 . ) • AMI * AM 1 *VV /2 . ) ) 
FF = (GAM-1. )*AMl*AMl/2. 

GG = l.-VV*VV 

AMD = VV<^AMl»SQRrF(l./(l.+FF*GG)) 

ZETAB = 1 .-VV*VV- (2./( ( GAM-1 . )«AM1*AM1 ))»(( (PP) **EXP1 )-l . ) 

WlY = W1*VV 

PlY = P1*PP 

TIY = TR1-C3*W1Y*W1 Y 

PRP » (TR.1/T1Y)*»EXP2 

PRIY = P1Y*PRP 

END 
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SUBROUTINE CURVE 



AA 

BB 

cc 

I 

J 

K 

M 

N 

NKAX 

ORD 

T1 

T2 

X 

XA 



FORTRAN NAMES. EQUIVALENT SYMBOLS AND DEFINITIONS 
\ 



a 

b 



Tl/2 



The terms of equation ORD = a + bx + cx 



Subscripts tor ordinate data points 



Counter to locate three nearest data points 
to input abscissa 

Array selector 

Number of data points in selected array 

Ordinate computed, loss coefficient 

value desired 

Data abscissa interval 

One^half of data interval Used to select the 
three nearest points 

Abscissa in terms of three nearest points 
Input abscissa 
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SUBROUTINE CURVE! N,XA1 , ORD, T I . NMAX , YY ) 

DIMENSION YYIUtSO ) , NM AX ( U ) , T I{ 4 ) tY( 50 ) 

COMMON Cl ,C2,C3,:4,GAM,EXP1 ,EXP2, EXP3,P0,T0,R 
XA = XAl 
NMAX=NMAX(N) ' 

TI=TI(N) 

DO 306 I = 1,NMAX 
506 Y( I)= YY(N, I) 

M = XA/TI 
X = TI*FLOATF(M) 

M = M+1 > 

IF (M-1) 304,300,301 

304 PRINT 212 
RETURN 

300 I = M 

J = M+1 
K = M+2 
GO TO 305 

301 NM = NMAX-1 

IF ( M-NM) 302,503,303 
303 I = M-1 
J = M 
K = M + 1 
GO TO 305 

302 DIFR = XA-X 
T2 = TI/2. 

IF (0IFR-T2) 303,303,300 

305 CONTINUE 
AA = Y( I) 

CC = 0.5« (Y(K)+Y( n-2.*Y{ J)) 

BB = Y( J)-Y(I )-C: 

XA = (XA/TI )-FLOATF( I-l ) 

ORD = AA+B8*XA+C: »XA*XA 

212 FORMAT (36H ABSCISSA WRONG IN SUBROUTINE CURVE ) 
END 
END 
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SUBROUTINE CPHI 



A 

CPHI 

DAB 

DAN 

DIFF 

EN 

EXPl 

EXP4 

EXP6 

EXP7 

GAM 

PHI 

PHI 2 

PHICR 

PPCR 

RC 

RA 

RB 

ZETA 

ZETAP 



FORTRAN NAMES. EQUIVALENT SYMBOLS AND DEFINITIONS 



n 






<!> 




(P/p) 



cr 



f 

fp 



Interim critical pressure ratio while computing 

Compute Subroutine name 

Name for the term (p/P)^^^^ 

Name for the term (p/P)^^^^ 

Difference between given and computed loss 
coefficient 

Polytropic exponent 

Exponent #1 = ( ^ - 1)/ ^ 

Exponent #4 » (1 - ^ ^ ^ 1)/^ 

Exponent #6 = 2/n 

Exponent #7 = (n+l)/n 

Specific heat ratio 

Flow function computed 

Flow function, squared 

Choked flow value of flow function 

Choked flow value of pressure ratio 

Computed loss coefficient 

Name for the term (p/P)^^^^ 

Name for the term (p/P)^^^^ 

Loss coefficient 

Polytropic loss coefficient 
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SUBROUTINE CPHI ( ZET A, A , PHI , PHICR , PPCR) 

COMMON Cl ,C2,C3,:4,GAM,EXP1 ,EXP2f EXP3,P0,T0,R 
ZETAP = -.001 
502 ZETAP = .OOUZETAP 

EXPU • ( 1 .-ZETAP) *(GAM-1 . )/GAM 
RA = 1 ./( A*»‘?XPli) 

= l./(A*»EXPl) , 

RC= ( RA-RB)/( 1 .-RB) 

DIFF = ZETA-RC 
IF (DIFF) 501,503,502 
501 ZETAP * ZETAP-. 03001 

EXPU = (l.-ZETAP)*(GAM-l.)/GAM V? 

RA * l./(A**EXP4) 

RB = l./(A**EXPl) 

RC= ( RA-RB)/( l.-RB) 

DIFF * ZETA-RC 
IF (DIFF) 501,503,500 
500 EN * GAM/(l.+ZETAP*(GAM-l.n 
EXP6 = 2./EN 
EXP7 = (EN+1. )/EM 
DAN = l./( A*«£XP6.) 

DAB = 1 ./ (A**EXP7 ) 

PHI2 » 2.*(GAM/(GAM-1 .) )«(DAN-OAB) 

PHI = S0RTF(PHI2) 

PHICR * ( (2./(EN+l.) )**( l./(EN-l.)) )*SQRTF(2.»GAM*(EN-1.)/ 
1 ((GAM-1. )*(EN+1.))) 

PPCR » ((EN+1.)/2.)**(EN/(EN-1.)) 

END 
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PROGRAMS TURPLOT AND TURBINE 



AKE 

AKIS 

ALFAl 

ALFA2 

BETAl 

BETA2 

DBETA 

DETA 

EFF 

ETA 

ETAA 

ETAC 

EXX 

NALFAl 

NPHIE 

PHI 

PHIR 

PSI 



FORTRAN NAMES, EQUIVALENT SYMBOLS AND DEFINITIONS 






Leaving loss coefficient 



Head coefficient, eq. 1 

Angle of flow exiting stator (radians), 
absolute 



0(2 Angle of flow leaving rotor (radians), 

absolute 






1 



Flow angle entering rotor (degrees), 
relative 



/2 

y 

7 

K 

% 

X 

^1 

f 

f 



Flow angle leaving rotor (degrees) , 
relative 

y ■ /^i • A 
'"7 - 7 • 7 a 

Reheat factor, defined in Ref. 1 

Efficiency 

Efficiency assumed 

Efficiency computed as an interim step for 
comparison 

Expression defined by eq. 5a 

Angle of flow exiting rotor (degrees) , 
absolute 

Carryover coefficient of kinetic energy into 
next stage 

Velocity coefficient in stator 

Carryover coefficient of kinetic energy from 
stator into rotor 

Velocity coefficient in rotor 
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PSIA 

RADRAT 

RBETAl 

RBETA2 

RELRAT 

RSTAR 

VELMl 

VELRAT 

WSQRl 

WSQR2 

XSQR 



• k 

A 









W ^ 
’*1 



w, 



Velocity coefficient in rotor computed for 
comparison 

Radius ratio through rotor 

Flow angle entering rotor (radians), relative 

Flow angle leaving rotor (radians), relative 

Ratio of relative velocities 

Degree of reaction, that fraction of the 
expansion through the stage which occurs in 
the rotor 

Meridional component of absolute velocity 
leaving stator/Peripheral velocity of rotor 
at mean diameter 

Velocity ratio of meridional components 
The relative velocity entering rotor, squared 
The relative velocity leaving rotor, squared 
Expression defined by eq. 5a 
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PROGRAM TURPLOT 

DO 50 NALFAl = 60 ,75,5 * 

DO 50 I a 1,12 V , 

RSTAR»-.05^(FLOATF( in*.05 
IF(N-2)21 ,22,22 • • 

21 PRINT kO 

N = N+1 ■ 

GO TO 23 

22 N = 1 . • 

PRINT 47 

23 CONTINUE 
NPHIE * 0 
RADRAT = 1.00 
VELRAT =1.00 
PHI = 0.95 
PHIR = 0.95 
EFF = O.C 

ETAC = .40 • ■ 

PS I * .99 

ALFAl = FLOATFINMFAl )*3. 1415927/180. 

AKIS » 2.0 
• J = 1 

13 AKIS = AKIS +2.0 
J = J + 1 . 

170RBETAl.a ATANF(TANF(ALFAT) - (( 1 . 0 )/( PHI »COSF ( ALFA 1) «SQRTF 
1 (1. - RSTAR)»S3RTF(AKI$) n ) 

BETAl * RBETAU130./3. 1415927 
12 CONTINUE 

OXSQR * (RSTAR*( 1. +EFF) ) ♦ (( PHI*PHI )*( 1 .-RSTAR 1 *( PHIR- ( VELRAT* 

1 VELRAT*C0SF(ALFA1 )*C0SF(ALFA1 J/(PSI*PSI n n - (2.*PHIR*PHI 

2 *SINF( ALFAl) *SQRTF( 1 .-RSTAR) /SQRTF ( AKI S ) ) + ( ( PHIR«>RAORAT* 

3 RAORAT-1. )/AKIS ) 

19 EXX = SQRTF ( XS3R ) . 

RBETA2*ATANF(-PSI *EXX/ ( PH I *VELRAT*COSF ( ALFAl ) *SORTF ( 1 .-RSTAR) ) ) 
BETA2 a RBETA2»130./3. 1415927 

20 DBETA a BETA1-BETA2 

PSIA a 0.99 - .000223*DBETA - 4. 97/ ( 1 80 .-OBETA ) 

DPSI a ABSFIPSI -PSIA) 

IF ( .0000001-DPSI ) 10,11,11 

10 PSI a PSIA 

GO TO 12 . 

11 PSI a PSIA 

OETA a (2./SQRTF(AKIS) )•( (PHI*SINF(ALFA1 )*SQRTF( l.-RSTAR) )- 
1 (RAORAT*RADRAT/SORTF(AKIS) )+(PSI*EXX*RADRAT)) 

OAKE aAKIS*( (PHI*3HI*( 1 RSTAR ) *COSF ( ALFA 1) *COSF ( ALFAl ) *V ELRA T* 

1 VELRAT) + ( (PSI*PSI )*XSQR)-(2.*PSI*EXX*RA0RAT/SQRTF(AKIS) ) 

2 +( (RAORAT*RAO^AT)/AKIS) ) 

VELMl a PHI*SORT= ( l.-RSTAR)*SQRTF(AKIS)*COSF( ALFAl ) 

0ALFA2 a ATANF( ( (RAORAT/SQRTFI AKIS) )-(PSI*EXX) )/(PHI*VELRAT* 

1 COSFIALFAl )*S3RTF( l.-RSTAR) )) . 

ALFA2 a ALFA2*13D./3. 1415927 

OWSQRl a 1.+(AKIS*PHI«PHI»(1.-RSTAR))-2.*SQRTF(AKIS)*PHI» 

1 SQRTF( 1. -RSTAR )*SINF( ALFAl ) 

0WSQR2 a (PSI*PSI»AKIS)*( ( (1 .+EFF)*RSTAR)+(PHIR*PHI*PHI* 

1 (1 .-RSTAR) )-(2 .*PHIR*PHI*SQRTF( l.-RSTAR )*S INF (ALFAl ) )/ 

2 ( SQRTF ( AKI S) )M PHI R+R ADRAT*RA0RAT-1. ) /AKIS) 

RELRAT a SQRTF(W5QR2)/SQRTF( WSQRl ) 

IF (J-6) 14,14,13 
15 PRINT 47 
J = 2 

140 PR I NT 41 ,NALFA1, NPHIE, RADRAT, VELRAT,PHI, PHIR, PSI ,R STAR, AKIS, 
1ETA,AKE,VELM1 ,REL RAT , BETAl , BETA2 , ALFA2 
IF (ETA-ETAC) 50,50,13 

400F0RMAT( /I 18H1 Al PE R2/R1 • VM2/VM1 PHI PHIR PSI R* 

1 KIS ETA KE VMl/Ul W2/W1 BETAl BETA2 ALFA2 /) 

410F0RMAT ( 13 , I 6, F7 . 1 , F9. 1 , F8. 2 , F7. 2 ,2F7.3, F8. 3 , F9.5 , 2F7 . 3, 

1 F8.3,2F8.2,F7.2 ) 

47 FORMAT (/) 

50 CONTINUE 
END 
END 
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/ (PHI4C0SF(4LFA1)«SQRTF 



PROGRAM TURBINE 
00 50 N» 1,4 
t4ALFA1 * 80 - (5* N) 

PRINT UO 

00 50 I «1,. 12 

RSTARs-.05+(FL0ATF( I) )*,05 
IF (1-7) 18,20,13 
20 PRINT 40 . ■ 

10 CONTINUE 

J « 1 • • 

K = 1 
NPHIE = 0 

ALFA! * FL0ATF(NALFA1)«3.1415927/180, 

RADRAT » 1.00 
VELRAT = 1.00 
PHI = 0.95 
PHIR = 0.95 
AKIS»8.0 
EFF » 0.0 
PSI = .99 
ETAA = 0.0 

170RBETA1 » ATANF(TANF( ALFA1 ) - (( 1.0) 

1 (1. - RSTAR)*S^RTF(AKIS)))) 

BETA! * RBETAU130./3. 1415927 

12 CONTINUE 

OXSQR * (RSTAR*( 1.+EFF) )♦( (PHI*PHI )*(1.-RSTAR)*(PHIR-(VELRAT» 

1 VELRAT*C0SF(AL=A1)«C0SF(ALFA1 )/(PSI*PSI ) ) )- ( 2. *PH IR»PH I 

2 »SINF(ALFA1)»SQRTF(1 .-RSTAR)/SQRTF(AKIS) )♦( ( PHIR*RADRAT» 

3 RAORAT-1. )/AKIS ) ) 

6XX » SQRTF(XSQR) 

RB£TA2=ATANF(-PSI *EXX/ ( PH I*VELRAT*COSF ( ALF A1 ) *SQRT F ( 1 .-RSTAR) ) ) 
BETA2 = RBETA2«130./3. 1415927 
DBETA * DETA1-BETA2 

PSIA = 0.99 - .0D0228*0BETA - 4.97/ C 1 80. -DBETA ) 

OPSl = ADSFIPSI -PSIA) 

IF ( .0000001-DPSI ) 10,11,11 

10 PSI = PSIA 

K » K+1 • • 

IF (K-1000) 12,12,45 

11 PSI » PSIA 

OETA a (2./SQRTF(AKIS) )•( (PHI»SINF(ALFA1 )*SQRTF(1.-RSTAR) )- 
1 (RADRAT«RAORAT/SORTF(AKIS) )-KPSI»EXX*RADRAT)) 

IF (ETA-ETAA) 13, 14, 14 

13 OETA = ABSF(ETA-ETAA) 

IF (.0000001-OETA) 15,16,16 

14 CONTINUE 
ETAA a ETA 
DKIS a (-.1)**J 
AKIS a AKIS+DKIS 
GO TO 17 

15 J » J+1 

IF (J-5) 14,14,45 
43 PRINT 44 
GO TO 50 

45 PRINT 46 . • 

GO TO 50 

16 CONTINUE 

OAKE »AKIS*( (PHI«>HI*( 1 ;-RSTAR ) -CDSF ( ALFAl ) *COSF ( ALFA 1 )* VELRAT* 

1 VELRAT )♦( ( PSI* PSI ) *XSQR)-( 2.»PSI*EXX*RA0RAT/SQRTF( AXIS ) ) 

2 -M (RAORAT*RADRAT)/AKIS ) ) 

VELM1 a PHI*SQRT= ( 1 .-RSTAR)*SQRTF( AKIS)*COSF( ALFAl ) 

0ALFA2 a ATANF( ( (RAORAT/SORTF( AXIS) )-( PSI*EXX) )/( PHI*VELRAT* 

1 COSF(ALFA1)»SORTF(1. -RSTAR)) ) 

ALFA2 a ALFA2*183./3. 1415927 

0WS3R2 a ( PSI*PSI*AXIS)*( { (1.+EFF)*RSTAR)+(PHIR*PHI*PHI* 

1 ( 1 .-RSTAR) )-(2 .*PHIR*PHI*SQRTF( 1 . -RSTAR ) *S INF ( ALFAl ) )/ 

2 (SORTF(AXIS) )«• (PHIR+RA0RAT*RA0RAT-1. )/AXIS| 

owsorL;^] i^a.*lcmtp(AKis )*phx* 

1 SQRTF( 1. -RSTAR )*SINF( ALFAl ) 

RELRAT a SQRTF(W5QR2)/SQRTF(WSQR1) 

OPRINT 4 1,NALFA1 , N PH IE , RADRAT , VELRAT.PHI , PHI R, PSI , RSTAR, AX I S , 

1 ETA, AXE, VELM1 , RELRAT, BETA 1 ,BETA2,ALFA2 
L a 1 
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19 CONTINUE 

ETAR * ETA-.005 

24 AKIS=AKIS+.l 
K=1 

ORBETAl = ATAMF(TANF(ALFA1 ) - (( 1.0) / ( PHI *COSF ( ALFA 1 ) #SQRTF 
1 (1. - RSTAR)*S3RTF( AKIS)) ) ) 

BETAl = RBETAl*130./3. 1415927 

22 CONTINUE 

OXSQR = (RSTAR»(l.+EFF))+( (PHI*PHI )*(1.-RSTAR)*(PHIR-(VELRAT* 

1 V£LRAT*COSF( ALFAl )*COSFtALFAl )/(PSI*PSI n )-(2.*PHlR«PHI . 

2 «SINF(ALFA1 )«SO^TF( l.-RSTAR)/SQRTF(AKIS) )♦( ( PHIR+RAORAT* 

3 RAORAT-1. l/AKIS ) ) 

EXX = SQRTF(XSQR) 

RBETA2*ATANF(-PSI •EXX/(PHI»VELRAT*C0SF(ALFA1 ) *SQRT F ( 1 .-RSTAR ) ) ) 
8ETA2 = RBETA2* 130. /3. 1415927 
DBETA * 8ETA1-BETA2 

PSIA » 0.99 - .0D0228*08ETA - 4.97/1 180. -OBETA) 

DPSI = ABSF(PSI -PSIA) 

IF (.0000001-DPSI.) 21,25,25 
21 PSI=?SIA 
K=X + 1 

IF (K-1000) 22,22,45 

25 PSI«PSIA 

OETA * (2./SQRTF(AKIS) )*( (PHI*SINF( ALFA! )*SQRTF(1. -RSTAR) )- 
1 (RADRAT*RADRAT/ SQRTFI AKIS) ) ♦ ( PS I*EXX*RAORAT ) ) 

IF (ETA-ETAR) 23,24,24 

23 CONTINUE 
M=1 

55 AKIS*AKIS-.00*1 
M=M + 1 

IF (1000-M) 51,35,35 

35 CONTINUE 

ORBETAl * ATANF(TANF( ALFAl ) - (( 1.0) / ( PHI »COSF (ALFA 1 ) *SQRTF 
1 (1. - RSTAR)*S3RTF( AKIS) ) ) ) 

BETAl » RBETAl*130./3. 141 5927 •* 

OXSQR * (RSTAR*( 1. 4-EFF) ) + ( (PHI*PHI)*(1 .-RSTAR)*(PHIR-(VELRAT* 

1 VELRAT«C0SF( ALFAl ) *COSF ( AL FA 1 ) / ( PS I *PS I ) ) )- (2. *PHIR*PHI 

2 •SINF(ALFA1 ) *SQ^TF( 1. -RSTAR) /SQRTF( AKIS) )♦( (PHIR»RAORAT* 

3 RAORAT-1 . )/AKIS ) ) 

EXX - SQRTF(XSOR) 

RRETA2*ATANF(-PSI «EXX/ ( PH I*VELRAT*COSF( ALFAl ) *SQRT F ( 1 .-RSTAR )) ) 
BETA2 = RB£TA2*130./3. 1415927 
O0ETA = BETA1-BEFA2 

PSIA = 0.99 - .030228*DBETA - 4.97/( 1 80.-DBETA) 

DPSI = ABSFIPSI -PSIA) 

IF (.0000001-DPSI) 36,38,38 

36 PSI=PSIA 
GO TO 35 

38 PSI=PSIA 

OETA = (2./SQRTF(AKIS) )*( (PHI*SINF(ALFA1 )«SQRTF( l.-RSTAR) )- 
1 (RADRAT*RADRAT/SQRTF( AKIS) ) ♦ ( PS I *£XX*RADRAT ) ) 

IF (ETA-ETAR) 55,26,26 

260AKE =AKIS«( (PHI*^HI*( 1 .-RSTAR ) *COSF ( ALFA 1 ) *COSF( AL FAl )«VELRAT* 

1 VELRAT)+( ( PSI*PSI) *XS0R)-(2.*PSI*EXX#RA0RAT/SQRTF( AKIS) ) 

2 +( (RADRAT*RADRAT)/AKIS) ) 

VELMl = PHI*SQRTF (l.-RSTAR)*SQRTF( AKIS)*C0SF(ALFA1 ) 

0ALFA2 = ATANF((UAORAT/SQRTF(AKIS))-(PSI»EXX))/(PHI*VELRAT* 

1 COSF( ALFAl )*SQRTF( 1. -RSTAR) ) ) . 

ALFA2 » ALFA2«18D./3. 1415927 

0WSQR2 = (PSI*PSI»AKIS)*( ((1.+EFF)*RSTAR)+(PHIR*PHI*PHI* 

1 (1.-RSTAR))-(2.*PHIR*PHI«SQRTF(1.-RSTAR)*SINF(ALFA1))/ . 

2 (SQRTF( AKIS) )«■ ( PHI R+R ADRAT*RA0RAT- 1 . ) /AKI S ) 

OWSQRl = 1.+(AKIS*PHI*PHI*(1.-RSTAR))-2.*SQRTF(AKIS)*PHI* 

1 SQRTFI l.-RSTA^ )*SINF( ALFAl ) 

RELRAT = SQRTF(WSQR2)/SQRTF(WS0R1 ) 

OPRINT 41,NALFA1 , N PHIE , RAORAT , VELRAT , PHI , PH IR, PSI , RSTAR, AK I S , 

1 ETA, AKE, VELMl , RELRAT, BETAl , BE TA2, ALFA2 
L»L + 1 

IFIL-3) 19,42,42 
42 PRINT 47 
GO TO 50 
51 PRINT 52 
GO TO 50 
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4OOFORMAT(/110H1 Al PE R2/R1 VM2/VM1 PHI PHIR PSI R* 

1 KIS ETA KE • VMl/Ul W2/W1 BETA1 8ETA2 ALFA2 /) 

410FORMAT ( / I 3, I 6, F7 . 1 , F9. 1 , F8 . 2 ,F7. 2, 2F7. 3 , F8.3 , F9,5 , 2F7. 3, 

1 F8.3,2F3.2,F7.2 ) 

44 FORMAT (/35H FIVE ITERATIONS OF AKIS COMPLETED) 

46 FORMAT (/34H 1003 ITERATIONS OF PSI COMPLETED ) 

52 FORMAT (35H 1000 ITERATIONS .OF AKIS COMPLETED ) 

47 FORMAT (///) 

40 FORMAT (IHI) 

50 CONTINUE 
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PROGRAM FANN03 



FORTRAN NAMES. EQUIVALENT SYMBOLS AND DEFINITIONS 



A 

DAB 

DAN 

EE 

EN 

EXPl 

EXP2 

EXP3 

EXP4 

GAM 

PHI 

PH12 

RA 

RB 

R 



P/p 



n 



<!> 



ZETA 

ZETAP 



’f 

/p 



Pressure ratio 

Name for the term (p/P)^^^^ 

Name for the term (p/P)^*^^ 

Reheat factor defined by ( ^ )/(l "^) 

Poly tropic exponent. See Appendix I 
Exponent #1 ■ (2/n) 

Exponent #2 » (n+l)/n 

Exponent #3 ■ (1 - ^p) ( ” 1)/^ 

Exponent #4 ■ ( 'Jf' - 
Specific heat ratio 
Flow function 

Flow function squared. Intermediate step In finding 

Name for the term (p/P)^*^^ 

Name for the term (p/P)^^^^ 

Computed loss coefficient to be used for 
determining p. R ■ (RA-RB)/(1-RB) 

Loss coefficient ■ 1 - 

Polytropic loss coefficient defined In Appendix I 
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11 

12 



13 



102 



TAPE 

TAPE 



PROGRAM FANN03 
DIMENSION PHI ( 11) 
DO 22 I » 125,143 
XI » I 

GAM * XI/100. 
lOOFORMAT (26H1FL0W 
19H PRESSURE ) 
3100F0RMAT (30H1 ZETA 
19H PRESSURE ) 
2100F3RMAT (32H1 REHE 
19H PRESSURE ) 
WRITE OUTPUT TAPE 
WRITE OUTPUT TAPE 
WRITE OUTPUT TAPE 
DO 11 M = 1,11 
ZETA(M) * -.025+* 
FORMAT (18H RATI 
WRITE OUTPUT TAPE 
WRITE OUTPUT 
WRITE OUTPUT 
A * 1.02 
CONTINUE 
00 in N « 1,11 
ZETA(N) » -.026+. 
ZETAP(N) » -.001 
ZETAP(N) * +.00U 
EXP3 * (l.-ZETAPI 
EXPU « (GAM-1. )/G 
RA » l./(A*«EXP3) 
RB * l./(A«*EXP4) 
R * (RA-RB)/( l.-R 
DIFF » ZETA(N)-R 
IF (OIFF) 101,100 
ZETAP(N) * -.0000 
EXP3 = (l.-ZETAP( 
EXP4 = (GAM-1. )/3 
RA « l./(A»*EXP3) 
RB « l./(A*»EXP4) 
R * ( RA-RB)/( 1.-^ 
DIFF a ZETA(N)-R 
IF (OIFF) 101,103 
EN a GAM/( l.*ZET<\ 
EXPl « 2./EN 
EXP2 * (EN+1. )/EN 
DAN « l./(A**EXPl 
DAB a l./(A»*EXP2 
PHI2 a 2.*(GAM/(G 
PHKN) a SQRTF(P^ 
£E(N) a (ZETAP(N) 
WRITE OUTPUT TAPE 
WRITE OUTPUT TAPE 
WRITE OUTPUT TAPE 
FORMAT (F6.2.13X, 
IF (A-1.99) 17,19 
A a A+.02 
GO TO 13 

IF (5.9-A) 22,20, 
A a A+.20 
GO TO 13 
CONTINUE ^ 

FORMAT (8FT2.5) 
END 
END 



101 



100 



111 



15 

17 

19 

20 
21 
22 

400 



,ZETA( 1 1 ),ZETAP( 1 1 ),EE( ^ 



FUNCTION FOR GAMMA > FS, 

POLYTROPIC FOR GAMMA a 

AT FACTOR *E* FOR GAMMA 

3, 10, GAM , 

4, 310, GAM 
5,210,GAM 

025*FL0ATF(M) 

0 ZETA a 11F9.3/) 
3,12rZE.TA 

4. 12, ZETA 

5. 12, ZETA . 



025«FL0ATF(N) 

ZETAP(N) 

N) )*(GAM-1. )/GAM 
AM 



Q) 

,102 

1+ZETAP( N) 

N) )*(GAM-1. )/GAM 
AM 



B) 

,100 

P(N)»(GAM-1. ) ) 



) 

) 

AM-1.))*(DAN-DAB) 

12 ) 

-ZETA(N) )/( l.-ZETA(N) ) 

3. 15, A, PHI 

4.15, A,ZETAP 

5.15, A,EE 
1 1F9.5) 

, 19 



20 



II ) 

2 // 

F5.2// 
a F5.2// 
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PROGRAM FANN04 



FORTRAN NAMES. EQUIVALENT SYMBOLS AND DEFINITIONS 
A Interim critical pressure ratio while computing 



CPHI 




Compute ^ , subroutine name 


DAB 




Name for the term (p/P)^*^^ 


DAN 




Name for the term (p/P)^*^^ 


DIFF 




Difference between given and computed loss 
coefficient 


EN 


n 


Polytropic exponent 


EXPl 




Exponent #1 ■ ( '^ - l)/'3(' 


EXP4 




Exponent #4 ■ (1 - " l)/'3^ 


EXP5 




Exponent #5 ■ 2/n 


EXP6 




Exponent #6 ■ (n+l)/n 


GAM 


'S' 


Specific heat ratio 


PH 


(P/P)„ 


Critical pressure ratio name for storing away 
and printing 


PHI 




Flow function assumed for comparison 


PHI 2 




Flow function, squared 


PHIC 




Flow function computed in subroutine based on 
and 5 


PHICR 




Choked flow value of flow function 


PPCR 


(p/p)„ 


Choked flow value of pressure ratio 


PR 


p/p 


Pressure ratio 


R 




Computed loss coefficient 


RA 




Name for the term (p/P)®*^^ 


RB 




Name for the term (p/P)®*^^ 
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z 




Poly tropic loss coefficient name for storing away 
and printing 


ZETA 


f 


Loss coefficient 


ZETAP 


^P 


Polytropic loss coefficient 
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Pf^OGRAM FANN04 

DIMENSION ZETAd n,PRm),PH(ll),Z(l 1) 

PRINT UOO 
DO 499 J = 1 ,n 

499 ZETA( J)*-.O25+.025*FLOATF(J) • 

PRINT 401»ZETA 

DO 490 I * 125»1'40 
XI = I 

GAM= XI/100. 

DO 498 K=1,n 
A = 1.8 
PHI=0.2 

497 CONTINUE 

CALL CPHI (GAM,Z£TA(KJ .A.PHICfPHICRfPPCRtZETAP) 

IF (PHIC-PHI) 494,493t495 

495 CONTINUE 
PHI * PHIC 
A = PPCR 
GO TO 497 

494 DIFF » ABSF(PHIC-PHI) 

IF (OIFF-. 0000001 ) 493.493,496 

496 CONTINUE 
PHI = PHIC 
A = PPCR 
GO TO 497 

4.93 CONTINUE 

PR(K) = PPCR 
PH(K) * PHICR 
ZIK) = ZETAP 

498 CONTINUE 

IF( 1-133) 492,491 ,492 

491 PRINT 400 
PRINT 401.ZETA 

492 PRINT 402, GAM 
PRINT 412, PH 
PRINT 403, PR 
PRINT 404,2 

4000FDRMAT (90Hl CHO<£D FLOW VALUES OF THE FLOW FUNCTION, PRESSURE RAl 
110 AND POLYTROPi: LOSS COEFFICIENT //) 

401 FORMAT (11X7HZETA = 11F9.3//) 

402 FORMAT ( 7H GAMMA* F4.2) 

412 FORMAT (18H PHI * 11F9.5/) 

403 FORMAT ( 1 8H (PO/P)C » 11F9.5/) 

404 FORMAT ( 1 8H ZETAP * 11F9.5//) 

490 CONTINUE 

END 

SUBROUTINE CPHI ( GAM, ZETA, A, PHIC, PHICR, PPCR, ZETAP) 

ZETAP « -.001 

500 ZETAP=0.001+Z6TA> 

EXP4= (1. -ZETAP )*( GAM- 1.) /GAM 

EXPl = (GAM-1 . )/SAM c 

RA = l./(A**EXP4) ^ 

RB = l./(A**EXPl) 

R = ( RA-RB)/( 1 .-RB) 

IF(ZETA-R) 501,502,500 

501 ZETAP*-. 00001+ZETAP 

£XP4=( 1. -ZETAP )*( GAM-1. )/GAM 
EXPl = (GAM-1. )/3AM 
RA = l./(A**EXP4) 

RB * l./(A»*EXPl) 

R = ( RA-RB)/( l.-RB) 

IF(ZETA-R) 501,502,502 

502 EN*GAM/( 1 .♦ZETAP* (GAM-1. ) ) 

EXP5 * 2./EN 

EXP6 = (EN+1 . )/EN 
DAN » !•/ ( A*»EXP5 ) 

DAB = 1 ./(A**£XP6) 

PHI2 = 2.*(GAM/(3AM-1. ) ) *(DAN-OAB) 

PHIC * SQRTF(PHI2) 

0PHICR = ( (2./(EN+l. ) )**( l./(EN-l. )) )*SORTF( (2./EXP1) •( (EN-1. )/ 

1 (EN+l.m 

PPCR * ((EN+l.)/2.)**(EN/(EN-l.n 
END 
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100 

101 



200 



102 

103 

104 



BB 

CC 

00 

IF 



105 



106 

300 



PROGRAM AFTER 
AME * 1.37 
NALFA *75 

ALFA = FL0ATF(NALFA)»3. 1415927/180. 

GAMMA « 1.37 

EXP = (GAMMA -1.) /GAMMA 

EE = (2. /((GAMMA-1. )*AME*AME1) 

PRINT 409 

PRINT 410iAMEtNALFA, GAMMA 

PRINT 40‘V 

00 50C L * 1.26 

PRAT = 1.02 - .02*FL0ATF(U 

OALFA = 0.0 

M * 0 

M = M*1 

CALFA = (-.1 )*»(M-1 l + OALFA 
OALFAR = 0ALFA*3. 1415927/180. 

COSRAT = COSF(ALFA)./CCSF( ALFA-OALFAR) 

AA » GAMMA*AME*AHE*C0SF(0ALFAR)-C0SRAT*(GAMMA-1. )«AME*AME/2. 

» -GAHMA*AME*AME-1 . 

= COSRAT*( l.+( ( gamma- 1 . )«AME«AME/2 . 1 ) 

* EB*3B-4.*AA*CC 
(DC) 101,200,200 

VELRATl = (-BB+SQRTF(00n/(2.*AA) 

OPRATl = (COSRAT»( 1 . + ( (GAMMA- 1 .) *AME*AME/2. )*( 1 . -VE LRAT 1 *VELR AT 1 ) ) ) 
1 /VELRATl 

2A = 1.-(VELRAT1*VELRAT1 )/(l.+EE-EE«(PRATl*«EXP) ) 

OIFF = ABSF(PRAT1 -PRAT) . 

IF (OIFF- .0001 1300,300.102 
IF (PRATl-PRAT) 103,300,101 
M » M+1 

OALFA * (-. 1 )** (M-1 l+OALFA 
OALFAR = 0ALFA*3. 1415927/180. 

COSRAT * COSF(ALFA)/COSF( ALFA-OALFAR) 

AA = GAMMA*AME*AME«C0SF(CALFAR)-C0SRAT*(GAMMA-1. )»AME«AME/2. 

= -GAMMA*AME*AME-1 . 

= COSRAT»( l.+( ( GAMMA-1. )*AME«AME/2.n 

* EB*BB-4. *AA»CC 
(00) 100,105,105 

VELRATl = (-B0+SQRTF(00) ) /(2. «AA) 

OPRATl = (COSRAT*( l. + ( ( GAMMA-1 .) *AME*AME/2. )• ( 1 .-VELRATl *VELR AT 1 ) )) 
1 /VELRATl 

ZA = ).-( V£LRAT1*VELRAT1 )/(l.+EE-EE«( PRAT1«*EXP) ) 

OIFF » ABSF(PRATI-PRAT) 

IF (DIFF-.OOOl) 300,300,106 
IF (PRAT-PRATl) 100,300,104 

PRINT 412,PRAT,PRAT1,VELRAT1 ,0ALFA,2A,AM0,AMA 
GO TO 
FORMAT 



BB 

CC 

00 

IF 



401 
4020F0RMAT 



1 



403 

404 

405 

406 

409 

410 

411 

412 
500 



1 



AME 

FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
MO 

FORMAT 
CONTINUE 
ENO 
ENO 



B2-4AC 

A 



500 
(20H 
(60H 
) 

(6F10.6/) 

(46H PRAT 
(/E20. 10) 
(F7.4,4F10.4) 
(25H1 MACH 
(F10.3,I5,F10.3 
(66H PRAT 

MA /) 
(7F10.4) 



IS NEGATIVE 
B 



/) 



PRATl VELRATl 



COSRAT 



OALFA 



GAMMA 



ZA / ) 



ALFA 

/) 



GAMMA 



/) 



PRATl VELRATl 



OALFA 



ZA 



- 156 - 






i f 



# il 



IQ'lg 



• > IJiigQLl M» « •« • tJ 

mmrn M »• • I . m •wMMi • J[£a<W>| 



fti 







I 






•- .» 



liripi ^ ^4)4 ••* I I 

l••■l••i|•‘ 4«p •••••• 



4 P»li 

• # 



•«ilM 





“v* 







Appendix IV- Critical Mach Number of Flow Undergoing Polytropic 
Expansion. 



The critical pressure ratio for a flow undergoing a polytropic 
expansion, as developed in Appendix I, is 




(1,13) 



To find the maximum velocity corresponding to this pressure 
ratio it must be recalled that, for a poly tropic process 







T ~ C 



constant 



where C is a constant 



(IV, 1) 



(IV, 2) 



Therefore 




(IV, 3) 



Using the T-s diagram of Fig. 24 as a reference, it can be seen that 



0ff^ 



and 



^ j€_r_ 



Z^J 



r 



/ 



J r-i 



w 



For the critical velocity then, using (1,13) 



\4 = 2-$^ JL ^ 

‘t-l 



/-) > 



'ml 

'n 






which reduces to 



=r 



-JLfjrsL) 

\ /Ti -h f J 



The acoustic velocity is 












which reduces to 






Therefore 



M ‘ 

'k/i 



1 / ' 






— / 

r-/ 



(IV, 4) 



(IV, 5) 



'■) 



(IV, 6) 



(IV, 7) 
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Also developed in Appendix I was the relation for the poly tropic 
exponent n, in terms of the poly tropic loss coefficient, as 



= 






• 5 ^ ('if-O +! 



Making this substitution, and collecting terms 



ym, — / 

T-f 



•J’ -•»'-/ 



and 



= 






For example, for 






1.37 and 
and for 






0.10, M 

cr 

0.20, M 
cr 



.931, 

.854. 



(1,7) 



(IV, 8) 



(IV, 9) 

I 
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